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ABSTRACT 



Supernova surveys can be used to study a variety of subjects such as: (i) cosmology through type la supernovae (SNe), (ii) star- 
formation rates through core-collapse SNe, and (iii) supernova properties and their connection to host galaxy characteristics. The 
Stockholm VIMOS Supernova Survey (SVISS) is a multi-band imaging survey aiming to detect supernovae at redshift ~0.5 and 
derive thermonuclear and core-collapse supernova rates at high redshift. In this paper we present the supernovae discovered in the 
survey along with light curves and a photometric classification into thermonuclear and core-collapse types. To detect the supernovae 
in the VLT/VIMOS multi-epoch images, we used difference imaging and a combination of automatic and manual source detection to 
minimise the number of spurious detections. Photometry for the found variable sources was obtained and careful simulations were 
made to estimate correct errors. The light curves were typed using a Bayesian probability method and Monte Carlo simulations were 
used to study misclassification. We detected 16 supernovae, nine of which had a core-collapse origin and seven had a thermonuclear 
origin. The estimated misclassification errors are quite small, in the order of 5%, but vary with both redshift and type. The mean 
redshift of the supernovae is 0.58. Additionally, we found a variable source with a very extended light curve that could possibly be a 
pair instability supernova. 

Key words, supernovae - general 



1. Introduction 

, During the last couple of decades supernovae have been shown 
to be powerful probes of both cosmology and star formation in 
high-redshift galaxies. A number of surveys with various sci- 
entific goals have been conducted and proposed. Most of them 
have in common that multiple imaging epochs are used to detect 
the supernovae. In some cases the surveys also contain spectro- 
scopic follow-up observations of the detected candidates. The 
spectroscopic information makes it possible to easily charac- 

. terise the detected supernovae and to measure redshifts, but this 
comes at the cost of added telescope time. Whole-sky supernova 
surve ys have alread y starte d, e.g., Pan-STARRSl (Young et al. 
l2008h and PTF (Rau et al. 2009), and several more are being 
planned, e.g., SkyMapper and LSST. With the large number of 
expected supernovae it will not be feasible to obtain spectra for 
all of them. Photometric techniques will thus be required to char- 
acterise the detected sources. 

Core-collapse supernovae (CC SNe) are the result of a mas- 
sive s tar ending its life in an energetic explosion (see, e.g.. lSmart3 
I2009h . Because the lifetimes of these massive stars are short. 



their presence signals that active star-formation is taking place in 
the host galaxy. The star-formation rate density for a cosmic vol- 
ume, defined by the redshift range and field size, can be derived 
from the core-collapse supernova rates by making an assump- 
tion on how many stars out of the current star forming popula- 
tion that explode. This method provides an independent tracer of 
the star-formation history o f the universe a n d has been used b y 
Dahlen et al, (2004j:,Botticella et al.l(l2008h : lBazin et alJ<2009l) . 
While Li et al.l ( 201 lal) also provide a local supernova rate, their 
big sample of SNe is also used to study the host galaxy proper- 
ties and connections to supernova subtype. Most of these surveys 
target low-redshift SNe (z < 0.3), but with deep e nough obser- 
vation s the technique also works at higher z (e.g. i Dahlen et al.l 
l2004l) . For this kind of projects it is important that the obser- 
vations are obtained as a "blind" survey to limit the selection 
effects. 

The use of thermonuclear (la) supernovae as standardis- 
able candles has been instrumental in measuring cosmologi- 
cal paramet ers for the ACDM concordance cosmology (e.g . 



logy (e^. 

aH I20T0I 

Wood-Vasev et al.li200'7HKessler et alJl2009 ). The systematic er- 



Astier et al.( [20 06; Ri ess et al.l [2 007; Amanu llah et 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile, under ESO 
programme ID 167.D-0492. 



rors resulting from the determination of the distance modulus 
for thermonuclear supernovae have been extensively studied and 
minimised. To achieve the precision needed to pinpoint cosmo- 
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logical parameters, it is necessary to obtain spectra for the su- 
pernovae. Despite the amount of work made on calibrating the 
distance modulus relation, the underlying physics of thermonu- 
clear supernovae is still partly unknown. In particular, the type 
of progenitor s ystem that g i ves ris e to the supernovae is uncer- 
tain; see, e.g., iRuiter et al.l (l2009l) for a discussion of the dif- 
ferent possibilities. One way of putting constraints on progen- 
itor models is to study the delay times of thermonuclear super- 
novae (i.e., the time between formation of the progenitor star and 
the supernova explosion). The rate of thermonuclear supernovae 
compared to either the global jP ahlen et al. 2004) or the local 
dSullivan et al.ll2006t iTotani et al 2008; Maoz & Badenesll2010l) 
star formation history can provide estimates on the delay time 
distribution even without observed spectra of the supernovae. 
A number of surveys targeting la SNe in galaxy clusters have 
also been undertaken ([ Sharon et al. 2010; Barbarv et al. 2010). 
Recently, iDildav et al.l (2010) presented an accurate measure- 
ment of the z ~ 0. 1 thermonuclear SN rate from the Sloan Digital 
Sky Survey II Supernova Survey. 

A major challenge to all intermediate- and high-redshift CC 
SN searches is how to take into account the effects of SNe m issed 
owing to large host g alaxy extinctions. iMattila et al.l (l2007l) and 
iKankare et alj ( l2008h present supernovae found in the nuclear 
regions of luminous infrared galaxies (LIRGs) by using infrared 
adaptive optics imaging. These types of searches are important 
to constrain the numbers of supernovae lost in LIRGs, in which 
most of t he massive stars are f ormed at intermediate and higher 
redshifts (iMagnelh et al.ll2009l) . 

Typing of supernovae using broad band colours in 
multiple epochs ha s been dem onstrated t o work by sev- 



eral a uthors (e.g., Barris & Tonrv. 20041 Johnson & CrottsI 



2006t iKuznetsova & Connolly! 120071: iPoznanski et alJ l2007d: 



SakoetalJ l2008r The application of the codes vary, in some 



cases the codes are used mainly to reject core-collapse SNe 
from the follow-up target lists of cosmological la surveys, but 
in others the codes are used as the only way of typing SNe 
without spectra. Most of the codes use template-fitting meth- 
ods where a number of k-corrected thermonuclear and core- 
collapse SN light curves are compared to the observed light 
curve and colour evolution. The simplest way of doing this fit- 
ting is through fitting. However, there are some problems 
with this approach because it does not allow prior information 
on, e.g., redshifts and peak luminosity to be used to full ex- 
tent. Kuznetsova & Connollv (2007) introduced a Bayesian ap- 
proach to supernova typing where probability distributions for 
the parameters could be used as priors . Lately, B ayesian methods 
have been further developed by Rodn ev & Tonr v (2009), where 
"fuzzy" templates are used to improve classification of SNe with 
non-standard light curves. With this method the templates are 
assigned an uncertainty that enters into the likelihood calcula- 
tions, which improves the classificati on quality of SNe with non- 
standard light curves. Recently, Kess ler et al.l (120 10a) presented 
the results of a supernova classification challenge, where a com- 
mon sample of supernova light curves was typed by different 
codes. Their results indicate that Bayesian typing codes are com- 
petitive and give reliable type determinations. 

The Stockholm VIMOS Supernova Survey (SVISS) is a 
multi-band (R + I) imaging survey aiming to detect supernovae 
at redshift ~0.5 and derive thermonuclear and core-collapse su- 
pernova rates. The supernova survey dat a were obtained ove r 
a six-month period with VIMOS/VLT (iLeFevre et all l2003h . 
[Melinder et al. (2008) describe the supernova search method 
along with extensive testing of the image subtraction, supernova 
detection and photometry. In this paper we report the discovery 



of 16 supernovae in one of the search fields of the survey and 
provide light curves along with type classifications for them. We 
will present the supernova rates and conclusions to be drawn 
from them in a future paper (Melinder et al, in prep.). 

The first part of the paper contains a description of the data 
set and the methods used to reduce it. In Section |3] we de- 
scribe the supernova typing method and the template supernova 
light curves. Section |4] contains a description of how the sim- 
ulations of supernovae have been set up and offers the results 
from extensive testing of the typing method. Finally, we present 
the results and possible implications in Section |5] and conclude 
by discussing and summarising the results in Section |6] The 
Vega magnitude system and a standard ACDM cosmology with 
[ho, Q.M, ^a} = {70, 0.3, 0.7) have been used throughout the pa- 
per. 



2. The data 

2.1. Observations 

The data were ob tained with the VIMOS instrument 
dLePevre et al.l l2003h mounted on the ESO Very Large 
Telescope (UT3) at several epochs during 2003-2006. The 
VIMOS instrument has four CCDs, each 2kx2.4k pixels with 
a pixel scale of 0.205"/pxl, covering a total area of 4x56 sq. 
arcmin. The SVISS observations were obtaine d in two fields, 
covering parts of the Chandra Dee p Field-South (iGiacconi et al.l 
1200 Ih and the ELAIS-Sl field dLa Franca et al.1 12004 . The 
observations in the ELAIS-Sl field were obtained in five broad 
band filters {U, B, V, R and /) centred at a ^ 00:32:13, 6 = 
-44:36:27 (J2000). The data used in this paper are only from 
the ELAIS-Sl observations, the CDF-S observations will be 
presented in a subsequent paper 

The supernova search filters are R and /, with roughly twice 
the exposure time in / compared to R. Observations in these 
filters have been divided into one reference epoch (henceforth 
epoch 0), seven search epochs and one control epoch. The im- 
ages used to construct the reference epoch were obtained in 
August 2003, the search epoch images were obtained from July 
2004 to January 2005, while the control epoch was taken in 
January 2006. The search epochs were separated by roughly one 
month. 

The UBV observations were obtained at several different 
epochs during the time period 2004-2006. Details on the data 
reduction and calibration of the UBV data along with a galaxy 
catalogue for the field will be presented in Mencia-Trinchant et 
al. (in prep). In the present paper the UBV observations have 
only been used to calculate photometric redshifts for supernova 
host galaxies. 

Table \T\ contains the specifications of the different epochs 
and parts of the data. It is worth noting that the quality of data 
is very good overall with median seeing of 0.76"and 0.68"in 
R/I, respectively, and all search epochs have seeing below 0.9". 
As can be seen in the table, the observations are also very deep, 
with mean 3cr limiting magnitudes of 26.8 and 26.3 in R and 
/, respectively. In Dahlenetal. (2008) the 50% detection effi- 
ciency magnitude of the survey in the F850LP filter is given, 
converting this to our magnitude system we obtain m/ ~ 26.5. 
This magnitude can be compared to our 3cr limits, based on the 
findings in Paper I; the depth is thus quite similar. The search 
for variable objects in the Subaru/X MM-Newton Deep Survey 
(iMorokuma et al.ll2008l; iTotani et al.l l2008) has a limiting mag- 
nitude of nij ~ 26.0, again converted to our magnitude system. 
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Table 1. Ovei-view of the ELAIS-Sl observations 



Filter 


Ep/Date 


Exp. time (s) 




Seeing (") 


R 


0/2003-08 


8 640 


— 


0.65 


R 


1/2004-07 


5 760 


26.83° 


0.68 


R 


2/2004-08 


8 280 


26.97° 


0.67 


R 


3/2004-09 


5 520 


26.71" 


0.69 


R 


4/2004-10 


8 760 


26.94" 


0.89 


R 


5/2004-11 


7 560 


26.80" 


0.80 


R 


6/2004-12 


7 320 


26.75" 


0.80 


R 


7/2005-01 


5 760 


26.88" 


0.76 


R 


C*/2006-jan. 


8 760 


- 


0.72 




0/2003-08 


11 520 




0.59 




1/2004-07 


1 1 520 


26.60° 


0.59 




2/2004-08 


18 240 


26.41" 


0.64 


J 


3/2004-09 


9960 


26.39° 


0.71 




Anr\f\A 1 n 


1 1 ^on 

i i JZU 


ZO. JZ 


U.Oo 




5/2004-11 


15 960 


25.85° 


0.76 




6/2004-12 


19 080 


26.04° 


0.64 




7/2005-01 


12480 


26.30° 


0.70 




C*/2006-jan. 


11520 




0.69 


U 


2003-2006 


37 740 


27.44^ 


0.73 


B 


2003-2006 


34 940 


28.53^ 


0.71 


V 


2003-2006 


20060 


26.93'' 


0.71 



Notes. The upper part of the table shows the observational properties of 
the Rl supernova search epochs, the lower part shows the properties of 
the UBV observations. The entries in the table are mean values over the 
four VIMOS quadrants. 

^°* 3cr magnitude limits estimated from the photometric scatter of sim- 
ulated sources placed in faint galaxies (m/ > 24.0) in the individual 
subtracted epoch images (see Section l4Tb . Control epochs, taken ap- 
proximately one year after the final SN survey epoch. 3cr magnitude 
limits using a 4xFWHM diameter aperture. 



Compared to other SN surveys, our data set is thus among the 
deepest ever obtained, although smaller in field size than others. 

2.2. Data reduction of the R and I observations 

The data were reduced using a data reduction pipeline writ- 
ten in MIDAS scripting language for SVISS developed by our 
team. Each image was bias-subtracted using an epoch master 
bias frame constructed by median-combining > 10 bias frames. 
Flatfield calibration data were obtained for each night in both 
filters and all science data were flatfielded using high signal-to- 
noise stacked frames. Cosmic-ray rejection routines were used to 
detect possible cosmic rays, but no automatic corrections were 
applied. Parts of images with suspected cosmic rays were man- 
ually inspected and pixels with cosmic ray contamination were 
flagged. Atmospheric extinction corrections were applied to the 
images using tabulated extinction coefficients from the ESO 
quality control web pages. The images were also normalised to 
counts per second. A bad pixel mask was produced for each of 
the images, containing vignetted (< 5% in most cases), saturated 
and cosmic ray flagged pixels. 

The VIMOS / band su ffers from quite severe effects of fring- 
ing (see'Berta et al.l l2008l for a detailed description of fringing 
for the VIMOS instrument), and care has to be taken to suc- 
cessfully remove these. We constructed a fringe map for each 
observation night by median-combining the non-aligned science 
frames and rejecting bright pixels by a standard sigma-rejection 
routine. To make this possible, the observations were made with 
a dithering scheme that was optimised to avoid having individ- 
ual science frames target the same area on sky. For some nights it 
was necessary to make two fringe maps owing to sky-brightness 



variations during the night. We did not need to use any mask- 
ing or smoothing on the resulting fringe maps, there were very 
few extremely bri ght source s in our field (causing optical ghosts, 
as mentioned bv Berta et al.r 2008) and the noise levels were in- 
significant compared to the individual frames. The fringe map 
was then subtracted from each of the flatfielded science frames. 
Finally, all of the science frames were manually inspected, some 
frames were removed because of vignetting and/or poor seeing 
(see Table[T]for the total exposure time used in each epoch). 

The good quality frames for each epoch were registered to 
a common pixel coordinate system using a shift- and rotate- 
transform (yielding a typical rms of ~ 0.1 pixels or 0.02"). 
The registered frames were then median combined. The bad 
pixel maps for each individual image were also registered and 
summed, yielding an exposure time map for each epoch. Some 
pixels (< 0.2% in the reference epoch) had an exposure time of 
zero, i.e., none of the individual frames contained useful data for 
that pixel. These pixels were flagged in the combined science 
frames. 

Each of the combined epoch images were photometrically 
calibrated using ~ 50 secondary photometric standards in the 
field. The calibrated photometry for these stars was obtained 
from stacked one-night images from nights where standard-star 
observations were available. The secondary standard sources 
were selected by requiring a signal-to-noise ratio of at least 100 
and a stellar-like point spread function (PSF). The sources were 
also required to be isolated in the single-night frame with no 
visible neighbours within a lOxFWHM (full width at half max- 
imum) radius. The photometry of the primary and secondary 
standards was obtained with IRAF/phot using an aperture size 
of lOxFWHM (with aperture correction). The resulting zero- 
point errors for each epoch are < 0.02 magnitudes, but it should 
be noted that all epoch images are scaled to the reference epoch 
zeropoint in the subtraction procedure. 

2.3. Supernova detections and pliotometry 

In' MeUnder et al.l (l2008 b. henceforth Paper I, we presented our 
supernova detection method, which is explained there in detail. 
We used IRAF/PYRAF scripts that were run in sequence and 
proceeded as follows: (i) accurate image alignment over the en- 
tire frame; (ii) convolving the better seeing image to the same 
PSF size and shape as the poorer seeing image, using a spa- 
tially varying kernel; (iii) subtracting the images; (iv) detection 
of sources in the subtracted frames, using both source detection 
software and naked-eye detection; (v) photometry and construc- 
tion of light curves of the detected transients. 

All search epochs were aligned to the pixel coordinate sys- 
tem of the reference epoch using the geomap/geotran tasks in 
IRAF/PyRAF. We found that the shift- and rotate-transform used 
when registering the individual frames was not good enough for 
registering the search images to the reference image. This is 
likely caused by changes in the geometry (e.g., differential re- 
fraction, changes in flexure of the telescope) of the frames over 
the time period our observations were obtained. To perform the 
registration at sub-pixel accuracy, we used a general transform 
that allows for shifts, rotations, shear, and pixel scale changes in 
the image being aligned. The reference sources for registration 
were bright, point-like objects in the field (approximately 100 
sources were used). The resulting standard deviations for the ge- 
ometrical transforms were smaller than 0. 1 pixels for all epochs. 

C onvolution t o a common PSF was made with the ISIS 2.2 
code dAlard & Lup ton 1998; Alard 2000). A convolution kernel 
was computed by comparing a number of reference sources in 
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the two images. The better seeing image (i.e., the image with 
smaller PSF width) was then convolved and scaled to match the 
photometry of the unconvolved image; the convolved frame was 
either the reference or the search frame. A background varia- 
tion between the frames was also computed and compensated 
for. The selection of suitable image sub traction par ameters for 
our dataset was investigated in Meli nder et al.l (1200 8). Finally, 
the reference epoch was subtracted from the search epoch and 
exposure time maps for the images were combined to make a 
weight map for the subtracted frame. 

2.3.1. Source detection 

Source detection was made in the subt racted frames using 
SExtractor (SE, Berlin & AmoutsI Il996b to obtain an initial 
source list. Separate detection was made in each epoch and fil- 
ter, except in the very last epoch (sources discovered in the last 
epoch would only have one point on the light curve and thus au- 
tomatically fail our selection criterion, see Section 12.3.21 ). The 
detection parameters for SE were set at a quite liberal level to 
make sure that no SN candidates were missed (i.e., this detec- 
tion threshold accepts sources at a lower signal-to-noise level 
than the rejection criteria, see Section I2.3.2t . Furthermore, the 
weight maps for the subtracted frames were used in source de- 
tection to lower the number of spurious detections. Using the 
weight maps reduced this number by about 20%, primarily close 
to the edges of the frames. 

The number of detections in a single subtracted frame is typ- 
ically in the order of 100 times the expected number of true 
varying sources. This is consistent with the findi ngs of other su- 
pemova surveys using similar techniques (e.g. Miknaitis et alj 
[2007). Most of these spurious detections are spurious subtraction 
residuals. About 50% of the spurious sources were then rejected, 
most of these by requiring that true sources must be present in 
both the R and / filter at the given epoch. Spurious sources close 
to saturated stars and image defects were also removed; the total 
number of candidates for all the epochs remaining after this ini- 
tial rejection procedure was ~1500. Most of them are spurious 
detections related to subtraction residuals of bright galaxies that 
are present in both filters. In Section|5]we show how constraints 
on the light curve and supernova typing can be used to safely 
reject the remaining spurious detections. 



2.3.2. Photometry 

Photometry on the detected sources was made using the IRAF 
daophot package. The PSF photometry was performed on all 
detected candidates using the task allstar. We also tried using 
aperture photometry with aperture corrections computed from 
the original worse seeing image, using the IRAF task phot. This 
does give fairly good results, but seems to be more susceptible 
to residual flux from the background galaxies than the PSF pho- 
tometry, thus giving larger errors (both statistical and systematic) 
in general. 

In Paper I we found that the photometric uncertainties esti- 
mated by daophot/phot were underestimating the true noise 
by a factor of two or more. The two main reasons for this are 
that (i) the pixels in the subtracted image will be positively cor- 
related owing to the convolution made in the image matching 
step; (ii) the sky noise is estimated in a region outside of the 
host galaxy, thus not taking subtraction residuals properly into 
account. Therefore, we used simulations to obtain reliable er- 
ror estimates. By simulating SNe at different brightness in each 



epoch and finding the scatter in their measured magnitudes, we 
obtained an estimate of the true photometric error for each point 
on the SN light curve. 

The simulations also allowed us to check the photometry for 
possible systematic errors. In Paper I we discuss the discovery 
of a small systematic flux offset (~I0% at the 3cr limiting mag- 
nitudes, lower for brighter sources) in subtracted frames. This 
offset is not present in all epochs/filters and changes between 
epochs. We can calculate this offset from the simulations and 
correct the photometry for this effect in the epochs where it is 
needed. 

R and / light curves are then put together for the 1459 sources 
remaining after automated detection in the subtracted frames. 
We rejected spurious detections by requiring the candidates to 
be brighter than the 3cr limiting magnitudes (see Table [T]) in (i) 
the detection epoch and the subsequent epoch and (ii) both the 
R and / filters. With this rejection criterion the number of super- 
nova candidates decreased to 115. At this point we are reason- 
ably certain that no true supernova candidates have been dropped 
(unless they were too faint to fulfil the rejection criterion), but 
there is likely still a number of spurious detections remaining. In 
some cases subtraction residuals will appear at the same location 
in both filters and in multiple epochs (e.g., very bright galaxies). 

2.4. Host galaxy redshifts 

Using the full set of UBVRI observations, we calculated pho- 
tometric redshifts for the supernova host galaxies. We used 
a template-fitting method with the mean galaxy lumi nosity 
function as a Bay esian prior, see Dahlen et al.l (I20I0 ) and 
iDahlen etall ( 120041) . We used 16 SEDs (spectral energy distri- 
butions) that were constructed by interpolating between four 
empirical templates (E, Sbc, Scd and Im galaxy types) from 
Coleman et al. (1980) and two starburst-galaxy templates from 
iKinnev et alj (II996 ). For the supernova host galaxies we used 
photometry from the stacked UBV images, but excluding frames 
obtained during the time period when the supernova is de- 
tectable. For the / and R images we used the reference epoch im- 
ages. These considerations ensure that the supernova light does 
not affect the determination of the photometric redshifts. 

In our observed subsection of the ELAIS-Sl field only two 
galaxies have known spectroscopic redshifts. This means that 
calibration and validating of the photometric redshifts is not 
straightforward. We obtained archive UBVRI data of the Hubble 
Deep Field-South (HDF-S) observed with the VIMOS instru- 
ment. The HDF-S has been observed extensively with spec- 
troscopy and we found 280 sources with spectroscopic redshifts 
ranging from z~0 to z ~3.5 that are also found in the VIMOS 
HDF- S observations (iVanzella et al.l 120021: iRigopoulou et al.l 
i2005t iGlazebrook et all l2006l) . We reduced and analysed the 
HDF-S photometric data in the same way as the ELAIS-Sl data 
set. The photometric calibration was also made in the same way 
(including corrections for galaxtic extinction) as for the ELAIS- 
Sl observations. 

Comparing the spectroscopic redshifts to the photometric 
redshifts obtained from the imaging data, we find a redshift un- 
certainty of Sz - 0.085 * (1 + z) and a frequency of catastrophic 
failures (defined as \Zphot - ZspeA > 0.3) of 9%. Because the in- 
strument and filters are the same, the analysis is made in exactly 
the same way and the depth of the imaging data is comparable, 
we assume that the uncertainties of the SVISS redshifts are the 
same as the uncertainties determined for this dataset. The un- 
certainty of the redshift determination was taken into account in 
the supernova typing, see Section 14711 The photometric redshift 
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technique and calibrations of it as applied to our data set will be 
presented in greater detail in Mencfa-Trinchant et al. (in prep). 

At first, the host galaxies were identified as the closest galax- 
ies to the SNe in terms of angular separation. But when the red- 
shift information was also considered, it became clear that some 
of the identified hosts were at unrealistical redshifts, either mak- 
ing any SNe in them too faint to be detectable or too bright com- 
pared to our templates. Instead of angular separation we there- 
fore used the photometric redshift information and chose the 
closest galaxy in terms of physical distance. With this method 
we were able to identify hosts for 13 of the SNe. For one of 
these (SN309) the chosen host was not the closest galaxy in an- 
gular distance. Identifying hosts for the SNe in this way is of 
course not unproblematic, there is a risk of selecting an incor- 
rect host galaxy. We did not try to estimate the percentage of 
misidentifications, but we reran the typing for all 13 SNe using a 
flat redshift prior and, in relevant cases, also with redshift priors 
based on other closeby galaxies. The result of these tests is that 
none of these SNe change main type and that the changes in red- 
shift are small. We conclude that host misidentification has only 
little effect on the main results of this paper 

For the remaining three supernovae (SN-14, SN-31 and SN- 
261) there are no nearby galaxies with redshifts that are consis- 
tent with hosting a SNe with the observed light curve. For these 
three SNe we ran the typing code without photometric redshifts 
supplied, thus using a flat prior on the redshift. This is also dis- 
cussed in more detail in Section 15.21 More information on the 
host galaxy identification and properties of these galaxies will 
be presented in an upcoming paper (Mencia Trinchant et al, in 
prep.). 

3. Supernova typing method 

Our typing method relies on a Bayesian template-fitting algo- 
rithm. The use of prior probabilities and Bayesian marginalisa- 
tion makes it possible to include previously known information 
on the different supernova types in the fitting technique while 
avoiding over-fitting. The goal of the method is to find the most 
likely supernova type (T) given an observed light curve {{F}) of 
a super nova candidate. The descr i ption of the Bayesian method 
follows iKuznetsova & Connolly! (l200 7h. although the notation 
and scope is slightly different. All calculations were made in a 
parallelised FORTRAN 90 code using double precision. 

Formally we want to find the type that maximises the prob- 
ability P{Tj\{Fi}), where - [1,...,A^7-] refers to the differ- 
ent types and / = [l,...,Ndp] to the points on the observed 
light curve. Nine supernova types (la, la/amt, ^^ibg-uke lagif-/*?, 
VoCnonnai, ^^Cbrighi, Hn, IIL and IIP) wcrc Considered. The to- 
tal number of data points {Ndi,) is given by the product of the 
number of filters {Nfnters) and the number of observation epochs 
(Nep). The SVISS observations were made in two filters {R and /) 
and in seven epochs, thus Ndp - 14 in this work. The probability 
P{Tj\{Fi]) cannot be calculated directly; using Bayes' theorem 
we can rewrite the probability: 



P{Tj\{Fi)) = 



P{{FiWj)P{Ti) 



(1) 



i:jP{\F,]\Tj)P{Tj) 

P{{Fi]\Tj) is the probability of obtaining the light curve data [Fi] 
for a given supernova type Tj. PiTj) contains all prior informa- 
tion on the supernova model light curve for the specific type. 
When Bayes' theorem is used in this fashion, we implicitly as- 
sume that the set of supernova types is complete, i.e., all super- 
nova candidates must be one of the considered types. For suf- 
ficiently peculiar supernovae, and non-supernovae. Equation [1] 



will not give valid probabilities. In Section l43] we describe how 
we use measures of goodness-of-fit to purge candidates with 
light curves that are too dissimilar from the model curves. 

For each of the nine supernova types we created template 
light curves using absolute magnitude {Mb) light curves and 
SEDs, which are described in more detail in Sections 13.31 and 
13.3.61 We used four parameters that uniquely define the light 
curve for a given type: (i) Mb the absolute rest-frame B band 
magnitude at peak, (ii) z, the redshift of the supernova, (iii) f, the 
time difference between the explosion date for the model and the 
first observational epoch,(iv) [Ry, E(B - V)] = ij; extinction in 
the host galaxy. We denote the template light curves by { fi j} 



{fij] = {fijiMB,z,t,n)], 



(2) 



and because these parameters uniquely define the light 
curve for a specific type, we may rewrite P{{Fi}\Tj) as 
Pi{Fi}\{fijiMB, z, t, /,))) and PiTj) as P(Mb, z, t, J/, Tj). 

We define the likelihood function for each supernova type: 

lj({FmB,z,t,n) = P{{Fi]\Tj)P{Tj) 

= P{{Fi)\[fi.j(MB, z, f, rimiMB, z, t, 1], Tj). 

(3) 

The probability that a supernova with a template light curve fij 
will have an observed light curve, {Fi}, is 



P({Fi}\{fij(MB,z,t,rj)})^Y] 



X f,-(Fi-fijf/(2SFf) 



^5Fi 



(4) 



where each observational data point is allowed to fluctuate ac- 
cording to Gaussian statistics (the widths of the distributions are 
given by the observational errors, 6Fi). It should be noted that 
the photometric quantities in this expression are in units of flux 
(counts/sec). Non-detections are included in the analysis as data 
points with zero flux and with an error given by the 1 cr limiting 
fluxes. 

The P{Mb, z, t, I}, Tj), or equivalently P(Mb, z, t, r}\T /), prior 
contains all prior information on the parameters for a given tem- 
plate. We assume that the parameters are independent and can 
thus factorise the prior: 



P{Mb,z, t, i]\Tj) = P(MB\Tj)P(z\Tj)Pit\Tj) 

xP{ij\Tj)P{Tj). 



(5) 



The individual parameter priors are described in the next section. 
The remaining prior, PiTj) contains information on whether a 
specific supernova subtype is more likely than others. The rela- 
tive rates of supernova subtypes are not well constrained at high 
redshifts. Therefore, we assume that the prior probability for a 
supernova candidate to be of a certain subtype is equal for all the 
types (i.e., a flat prior). 



Nr 



(6) 



There are more priors/parameters that could have been used 
to construct the template prior Including more prior informa- 
tion will of course affect the typing and can make it more accu- 
rate. For this to work, the parameter in question needs to have a 
known probability distribution (or at least a valid range of val- 
ues). It should be noted that adding more priors will cause the 
computation time of the typing to increase by a factor equal 
to the number of steps used for the probability distribution. 
Our choice of four priors is thus based on a compromise be- 
tween limiting computing time and choosing parameters with 
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well-known distribution on one hand and typing accuracy on 
the other. Other authors have used additional or a different set 
of priors. Kuznetsova & Connolly (2007) include prior informa- 
tion on stretch in their typing code, using a Gaussian distribution 
with mean and width from observations. In our code we use ther- 
monuclear supernova templates with different stretch to account 
for this variance. The inclusion of a colour uncertainty to the 
templates can also be desc ribed in terms of an additional prior, 
this has been used by, e.g.. iPoznanski et aTl (|2007a). We do not 
have colour uncertainty in our typing, instead we choose to use 
more templates that have different colour evolution. 



supernova type, thus we may write P(z\T) - P{z) and obtain 



g-(z-<>)V(2&') 

Piz) = — Az. 



(8) 



Az is the numerical step size, similar to AM in eq. Q The flat 
redshift prior is given by 



P{z) 



Az 

Zu - Zl ' 



(9) 



3.1. The parameter prior distributions 

Each of the four parameters has an associated probability distri- 
bution function. The parameter priors on the right hand side of 
equation |5] are given by the respective probabilities. 

The prior on the peak absolute B magnitude is given by 



P(.Mb\T) 



-AMr 



(7) 



where (Mb) is the mean peak magnitude for a given super- 
nova type (see Table |2] and the discussion in Section 13.31 1. ctm 
is the dispersion of the mean magnitude and AM is the nu- 
merical step size. Assuming the peak magnitudes to be nor- 
mally distrib uted around the mea n is consistent with observa - 
tions, c.f. Richardson et al. (2002) and Richardson et al. (2006). 
The range in peak absolute magnitudes considered is given by 
[{Mb) - 2a-M, (Mb) + 2o-m] for each subtype and the numerical 
step size used is ~ 0.02 - 0.04. 

The redshift prior used is either a Gaussian or a flat distri- 
bution, depending on whether the host galaxy has a valid pho- 
tometric redshift available or not. The mean, Zp and width, 6z 
of the Gaussian distribution is given by the photometric redshift 
fitting of the host galaxy (see sec. 12.41 and Mencia-Trinchant et 
al., in prep.). The redshift of the supernova is independent of the 



Table 2. Properties of the supernova photometric templates 



Type 


(Ms) 


ctm 


Stretch 


Fraction 


References 


Ia-91T 


-19.64 


0.30 


1.04 


0.125 


1-3 


la - normal 


-19.34 


0.50 


1.00 


0.615 


1,3-4 


la - faint 


-18.96 


0.50 


0.80 


0.196 


1,3-6 


Ia-91bg 


-17.84 


0.50 


0.49 


0.064 


1,3-4 


Ibc - bright 


-19.34 


0.46 


N/A 


0.057 


1,7-9 


Ibc - noiTnal 


-17.03 


0.49 


N/A 


0.259 


1,7-9 


IIL 


-17.23 


0.38 


N/A 


0.028 


1, 10-12 


Iln 


-18.82 


0.92 


N/A 


0.248 


1, 12-13 


IIP 


-16.66 


1.12 


N/A 


0.407 


1, 12, 14-21 



Notes. (Mb) is the absolute magnitude in the lohnson-B filter at peak 
(assuming standard ACDM cosmology {ho,Q.M,^A] = {70,0.3,0.7)), 
(Tm is the dispersion in the peak magnitude. Fraction refers to the ex- 
pected fraction in a magnitude limited sample. 



Refer ences. (l )lNugentl (1200' 
( 120031); (4) iNugent et al " 



( I2OO2E 



JhaetalJ (l2006h: (1) IL evan et al 
(9) ^Richardson et al.'' (I200Q) : 



): (2)IStem et al.l (l2004h: (3)lTonrv et al.l 
(im^; (5) Thillips et al.' (11999)) 



(in ICap pellaro et al. (199% 
('13)lDiCarloe t al. (2002) 
( 120041) : (16) Elmhamdi et al. (200 



(2005); (8) Hamuvetal 
(10) Gilliland et al. (1999); 
(12) Richardson et al. (200%; 
(14) Des sM-t et al.l ( 120081); (15 ) Baron et af 



(18)'Hendrv et al. (2005); (19) Sahu ltai] ( [2006h ; (20) lPastorello et al.l 
(I2OO6) ; f21) ,0uimbv et al. t2007.) 



(17) [ Leonard et al. (2002) 



where [zhZu] is the lower and upper limits of the redshift in- 
terval considered. For both distributions we used a numerical 
step size of 0.02 in redshift and a redshift interval given by 
[zi,Zu] = [0.0,2.0]. 

The probability of a given time difference is assumed to be 
equal for all time differences between f] to f2. The time differ- 
ence probabilities do no depend on the type, P(f|7~j) = P{t), and 
the prior is 

At 1 

^^(0 = = (10) 

ti - h N, 

where Af is the numerical step size and A^, the number of steps. 
We used Af = 1 day throughout this work. The time difference 
interval considered is given by [fi, fa] = [-136, 142], this inter- 
val is chosen based on the time difference between our search 
epochs and that we require supernova candidates to be observed 
in at least two subsequent epochs. 

The extinction m odels considered are two Cardelli laws 
(ICardelli et al.H 19891) with Ry = 2.1 and 3.1 and the Calzetti 
law (Calzetti et al. 2000) with Ry - 4.05. Some authors have re- 
ported indications of steep extinction laws for supernovae (e.g. 
Goobar 2008), this motivates the inclusion of the Ry - 2.1 law. 
The Calzetti law is normally used in star-forming galaxies and is 
therefore a reasonable assumption for core-collapse supernovae. 
However, all extinction laws are considered for any given type. 
The second parameter for the extinction is the colour excess, 
E{B - y), the range of excess considered is 0.0 < E{B - V) < 0.6 
with step size of 0. 1 . The total number of combinations of the 
different extinction laws and values of E{B - V) (r\) is Next -21. 
We assume a flat prior that does not depend on the type under 
consideration and arrive at 



P(il\Tj) 



1 

Next' 



(11) 



We performed tests on simulated and real data that indicate 
that the typing is fairly insensitive to changes of the extinction 
prior range. Using a higher maximum extinction or a different 
step size does not change the type of any of the observed SNe 
(see Section[5]). 



3.2. Determining tlie most likely supernova type 

When trying to find the most likely supernova type, the four pa- 
rameters are all nuisance parameters, we marginalize the likeli- 
hood function over the four parameters to get the summed like- 
lihood function, or in Bayesian terms, the evidence: 

= Z Z Z J]P({Fi}\{fij(MB,z,t,ii)mMB,z,t,i],Tj). 

Mb z I I] 

(12) 
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Combining this, we anive at the following formula for the evi- 
dence (given that photometric redshift for the host is available): 



A^^ Zl A?, Zl 1^^^^ Z 



-(Ms-<Mb))2/(2o-=„) 
-(Fi-f,jfK25Fj) 



-AMb 



(13) 



The probability P{Tj\{Fi]) is then equal to the relative evi- 
dence for each supernova type, 



P{Tj\{F,]) = 



(14) 



The most likely type for a given SN candidate is the one with 
the highest P(Tj\{Fi}). We co-add the probabilities for the sub- 
types belonging to either of the two main types (thermonuclear 
and core collapse types) and obtain the probabilities P(TN) and 
P{CC) for each typed supernovae. Our tests of the typing using 
simulated data of the same quality and type (see Section l4Tb 
indicate that the resulting subtype classifications within these 
two main types are prone to quite large errors (in the order of 
10-20%), while the errors on the main type classifications are 
significantly smaller (in the order of 5-10%). 

In Section |4] we discuss how to set constraints on the ev- 
idence and normalised probability to reliably reject both mis- 
classified and non-supemova objects. 

We also constructed a "best-fit" light-curve for each super- 
nova. This was made by fixing the supernova type and then find- 
ing the most likely value for each parameter one at a time, by 
marginalising over the remainin g three parameters. As noted by 
iKuznetsova & Connolly! ( |2007|) . Bayesian techniques do not al- 
ways give the best estimate for the individual fitted parameters. 
To obtain trustworthy estimates of individual parameters, dedi- 
cated simulations have to be run, which is beyond the scope of 
this paper We merely used the fitted light curves as a measure 
of the overall quality of fit (a method that is studied extensively 
by Monte-Carlo-simulated light curves, see Section l4Tt . 

3.3. Supernova subtype templates 

We used supernova spectra and absolute magnitude light curves 
to construct light curve templates for nine different supernova 
subtypes. The spectra and Mb light curves fo r all the s e type s 
were (with a few exceptions) obtained from iNugenfl (|2007|) . 
The spectra we used extend from 1 000 Ato 25 000 Aand cover 
epochs from day ~l-2 to ~100 after explosion. The light curves 
were scaled to the a bsolute peak magnitu des given in Table|2]and 
were obtained from lDahlen et al.l(l2004l) (and references therein) 
unless otherwise noted. The table contains a complete list of the 
references used in the construction of each template. 

The output light curve template library consists of observer- 
frame light curves for the redshift range 0.01-2.0 (with a step 
size of 0.01 in z) in the VIMOS R and / filters for each of the 
nine supernova subtypes. As an example. Figures [T] and |2] show 
the R band light curve and the R-I colour evolution for the nine 
templates at z = 0.5 (roughly corresponding to rest frame B and 
B - V). It should be noted that the similarities in colour between 
the Ibc-bright and Ibc-normal templates and for the IIP and IIL 
templates, respectively, arise because of the spectral templates 
used (see relevant sections below). 

The subtype f ractions li sted in Table |2] are based on the 
measurements of iLi et al.l (1201 Ibl) and the compilation by 



__ la-91bg 

Ia-91T 

— la-normal 

■ ■ ■ la-faint 




Epoch (days) 

Fig. 1. Template light curves for the four thermonuclear super- 
nova types at z = 0.5. The upper panel shows the R light curve 
for the la-normal subtype as solid (black), the la-faint subtype 
as dash-dotted (blue), the 91T-like subtype as dotted (green) and 
the 91bg-like subtype as dashed (cyan). In the lower panel the 
R-I colour evolution for the four subtypes is shown using the 
same line styles (this figure is available in colour in the electronic 
version of the article). 



iDahlen et al. I (l200l and were calculated for our magnitude lim- 
ited sample using the limiting magnitudes from the light curve 
rejection step. These fractions were not used in the actual typ- 
ing, but were only used to study the misclassification ratios, see 
Section 1431 



3.3.1. Type la supernovae 

We used four different light curve templates for the thermonu- 
clear supernovae. The four templates (91bg-like, la-faint, la- 
normal, 91T-like, from fainter to brighter) have different abso- 
lute peak magnitudes and consequently different decline rates 
as well. The decline rates were p arameterised usir i g the stretch 
parameter, s, first introduced by iPerlmutter etaP (11997). The 
range in stretch (s - 0.49-1.04) for the four templates provides 
a sparsely sampled grid that cov ers the observed stret ch distribu- 
tion of type la supernovae (e.g. .' Sullivan et af]|2006l) . 

The spectra and rest-frame B light curve used to construct 
the la-normal and la-faint templates were first presented in 
iNugent et al.l (l2002l) . The absolute peak B magnitu de used for 
the la- normal template, -19.34, was obtained from Tonrv et al.l 
(l2003h . The peak magnitude for the la-faint subtype was calcu- 
lated by starting from the la-normal peak magnitude and scal- 
ing to a 5 = 0.8 peak m agnitude using the luminosity-to-decline 
rate relationship from Phillips et al.l (1 19991) and the decline-rate- 
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Ibc-brighl 
Ibc-normal 

- ML 

- Iln 

- IIP 




Epoch (days) 



Fig. 2. Template light curves for the five core-collapse supernova 
types at z = 0.5. The upper panel shows the R light curve for the 
IIP subtype as thick solid (black), the Iln subtype as dash-dotted 
(red), the IIL subtype as solid (green), the Ibc-normal subtype 
as dotted (black) and the Ibc-bright subtype as dashed (blue). In 
the lower panel the R - I colour evolution for the five subtypes 
are shown using the same line styles (this figure is available in 
colour in the electronic version of the article). 



to-stretch relation from lJha et al.l (l2006b . The resulting absolute 
peak B magnitude is -18.96. 

The light curv e and spectr a for th e 91T-like template were 
first presented in IStem et all (l2004b . These input data were 
corrected for extinction assuming {Rv,E{B - V)} - {3.1/0.2) 
(Nugent 2007). The 91bg-like template is based on light curves 
and spectra of SN 1991bg and SN 1999by dNugent et al.ll2002l) . 
The pea k absolute B mag nitudes are, respectively, -19.64 and 
-17.84 ( Tonrv et al.l2003h . The templates have stretch values ac- 
cording to their namesakes, 1.04 for the 91T-like template and 
0.49 for the 91bg-Uke template. 



3.3.2. Type Ib/c supernovae 

We used two Ib/c supernova templates, the Ibc-bright and Ibc- 
normal templates. The spectral evolution is based on the work of 
iLevan et al l (l2005h and the light curve is based on observations 
of SN 1999ex from lHamuv et al.l(l2002h . The observed data were 
corrected for extinction because the used SN h ave notab le ex- 
tinction, {Rv,E(B- V)\ = {3.1,0.4) is assumed dNueenllbOOTl) . 
iRichardson et al.l (l2006h find that the absolute peak magnitude 
distribution of Ib/c SNe can be described by using one bright 
and one normal population. We used the mean absolute peak 
magnitudes, -19.34 and -17.03, and scatter for the two popula- 
tions from their paper; converting the magnitudes to our chosen 
cosmology and k-correcting to the B band. The data available 



in the SVISS do not allow a refinement of the Ibc typing into 
lb and Ic subtypes, and it is not clear-cut that these two spec- 
tral types correspond to the two different ph otometric types (the 
bright and norm al) considered in this work (iHamuv et al.ll2002l : 
iLi et al.ll20TTbh . 

It should be noted that the red colour of the template in pre- 
peak epochs (the first 5-7 days in rest-frame) is the result of bas- 
ing the t emplate on SN 199 9ex, which showed this effect at early 
epochs dHamuv et al.l2002h . There are not many Ib/c supernovae 
that have been observed before peak brightness, but observations 
of other Ib/c's (e.g., SN 2008D from Modjaz et al. 2009) show 
that this effect may not be characteristic of the type. We per- 
formed tests using the typing code with a Ib/c template without 
the early red phase and found that the choice of template has 
very little influence on the resulting probabilities (less than 1%). 



3.3.3. Type IIL supernovae 

The II L spectra and light curves were first used inlGilliland et^ 



1999 ). The light curve is or iginally from Cappellaro et aL 



1997 ). IRichardson et al.l (l2002h presented the absolute B peak 



magnitude distribution of IIL SNe, and we adopt the mean 
peak magnitude and scatter for their normal IIL population. 
Converting the magnitude to our cosmology, we arrive at a B 
peak magnitude of -17.23. 



3.3.4. Type Iln supernovae 

The spectra and li ght curves are based on SN 1999el 
dPi Carlo et al.ll2002h . whereas the absolute B peak magnitude 
and scatter are obtained from IRichardson eTal . (2002). After 
compensating for the different cosmological parameters, the re- 
sulting peak magnitude is -18.82. 



3.3.5. Type IIP supernovae 

The spectra used to build the IIP supernova template come 
from different sources. For the early epochs (< 33 days 
after exp losion) we us e d ext inction-corrected spectral mod- 
els from iDessart et al.l (l2008l) . which in turn are based on 
SWIFT UV-optical observations of SN 2005cs and SN 2006bp 
(Pastorello et al. 2006; Brown et al. 2007; Immleretal. 200^ 
lOuimbv et al.l I20p7j). The late epo ch spectra are ba sed on 
SN 19 99em from lBaron et al.l (l2004l) and obtained from .Nugent! 
(120071) . The reason for not using the SN 1999em spectra for 
the early epochs is that the UV part of those spectra are mod- 
elled and extrapolated from optical observations, the more recent 
modelling based on early epoch UV/optical observations should 
therefore provide a more accurate spectral evolution. 

The light curves were constructed using photometric data 
for SNe 1999em 1999gi, 2 003gd, 2004et 2005cs and 2006bp 
E003; Leo nard et al.ll2002l; iHeiidrv et al 20051; 



('Elm hamdi et al. - ^ , 

S ahu et al.l 120061; IPastorell o et al.1 l2006t lOuimbv et al.1 l2007h 
The light curve of each of these SNe was scaled to a common 
peak magnitude and corrected for extinction using the values 
given in the original references. Each curve was then re-sampled 
with a resolution of one day over the range of observed epochs 
(using spline interpolation). The final IIP light curve was then 
obtained by averaging over the six interpolated light curves. The 
adop ted absolute B pe a k mag nitude, -16.66, and scatter come 
from IRichardson et al.l (l2002i) . again compensating for the cos- 
mological parameter difference. 
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3.3.6. K-corrections of SN light curve templates 

We then calculate VIMOS R and / apparent light curves for 
the redshift range [0, 2] using the following formula (similar to 
[Pahlen e t al. 2004): 



my(t, z) ^Mpeak,B + ^Meit X (1 + z) ') + niz) 
+ AB-K^(z,tx(l+zr'), 



(15) 



where y refers to the observed filter (R or /), t is the observer 
frame epoch, Mpeak,B is the peak absolute magnitude in a rest- 
frame B filter, AMg is the light curve decline relative to the peak, 
fi{z) is the distance modulus and is the rest frame extinction in 
th e host g alaxy. The K-correction K^, following the formalism 
of iKim et al. ( 1996), is given by 



K^(z,t)^- 2.5 log 



+ 2.5 log 



+ 2.5 log(l +z) 



/ Z(A)SB(A)dA 
J Z(A)Sy(A)dAj 

/ F(A,T)SB(A)dA 
jF{A/{l+z),T)Sy(A)dAj 



(16) 



In this expression Z(A) is the spectral energy distribution of 
Vega; S b(A), S ,j(A) the Johnson B and VIMOS R/I total filter 
transmission curves and F{A, r) the spectral energy distribution 
of the supernova at rest frame epoch r. 

4. Testing the typing method 

4.1. Simulated supernova light curves 

To test the typing accuracy of the code applied to the observed 
SNe, we created a simulated mock detection catalogue. This cat- 
alogue was created to resemble the expected characteristics of 
the detected SN sample when it comes to peak magnitude distri- 
butions, light curve shapes, redshift distribution, and dust extinc- 
tion. The catalogue also reflects the specifics of the SN search, 
including the cadence between observing epochs and the limit- 
ing magnitudes in the detection filters. 

We first assumed a total SN rate by summing th e ther- 
monuc lear (TN) and core-collapse (CC) rates from Dahl en et al.l 
(12004 . We used this rate to assign a random redshift to each 
mock SN in the range < z < 2. The SNe were thereafter given 
a main type, either TN or CC, from the relative strength of the 
rate of these types at the assigned redshift. The thermonuclear 
SNe were then subdivided into a faint, a normal, a 91bg-like and 
a 9 It-like population, while the CC SNe were subdivided into 
IIL, IIP, Iln, Ibc-bright or Ibc-normal. It should be noted that 
the fractions of each simulated subtype are equal at this point. 
We did this to ensure that the number of simulated SNe for each 
subtype is sufficiently high for statistics, even for the rare sub- 
types. The fractions given for each subtype in Table|2]were used 
when analysing the typing results for the simulated light curves 
(see Section l43] l. 

A peak B-band magnitude was assigned to each SN 
(Table |2ll. This was then perturbed using the peak magnitude 
dispersion values given in the table. Each SN was assigned an 
explosion date over the course of one year We then used the ac- 
tual cadence of the SN search and the B-band light curves for 
each specific type to calculate the absolute magnitude of the SN 
at the different observational epochs. Because the SN explosion 
dates are distributed over time, we will have some SNe that are 
observable in all epochs (i.e., those exploded before the first ob- 
servation), and others that explode during the observational pe- 
riod and thus only observable in some of the epochs. 



In the next step we used the type specific SEDs of the SNe to 
calculate the K-corrections that give us the apparent magnitudes 
in the observed R and I filters corresponding to the absolute B- 
band magnitudes (see Section [3.3.61 ). The apparent magnitudes 
were also corrected for extinction in the SN host galaxi es using 
the extinction distributions from iRiello & PatatI (l2005l) . In the 
final step, we added a photometric error to the apparent magni- 
tudes. This consists of both a statistical part, derived from the 
expected S/N=3 limits in R and I (see Section l2.3.2l and a "sys- 
tematic" part, for which we assumed an extra 4% error in the 
flux. The simulated catalogue consists of apparent magnitudes 
and errors in R and I at each observational epoch. We treated 
objects with S/N< 1 as non-detections. This is also the case for 
epochs observed prior to the explosion of a particular SN. 

Because we rely on photometric redshifts, we also assigned a 
simulated redshift to each SN, Zsim- We calculated this by adding 
a random scatter, Azmd - drawn from a Gaussian distribution 
with CT; - 0.06 - to the true redshift (z,rue) according to Zsim = 
Ztrue + (1 + Zrme) * Azrnd- The (T, uscd here is the photometric 
redshift accuracy of our host galaxy catalogue, see Section 124] 

When testing our code using the mock catalogue, we culled 
the sample by applying additional selection criteria. As our de- 
fault set-up, we only included objects with a S/N>3 in at least 
two consecutive epochs in both R and I (mimicking the set-up 
used for the observed supernovae). The total number of remain- 
ing simulated supernovae after the cull is ~ 18 000, which corre- 
sponds to roughly 2 000 SNe per subtype (see Table|3). 

4.2. Rejection of non-supernovae 



As mentioned in Section 13.21 the Bayesian typing method only 
gives valid results when the input light curve can be modelled by 
one of the templates. If the variability of a source has a non-SN 
origin (or, though less likely, is a SN with a very irregular light 
curve), the resulting likelihoods will be very small. Indeed, the 
evidence will be zero (within machine accuracy) in many cases 
where the observed light curve is just too different from a tem- 
plate SN light curve. The evidence, (see Equation [TlTl, can 
therefore be used to reject possible anomalous sources (AGN, 
variable stars, co mplicated subtraction resi d uals). We followed 
the suggestion of iKuznetsova & Connolly! (l2007h and defined 
rejection limits for each template, X99.9,r- These were calcu- 
lated by finding the evidence below which fall 99.9% of the 
evidences - resulting from fitting of the Monte-Carlo-simulated 
light curves. The likelihood thresholds are in the order of 1 

It should be noted that the likelihood threshold is only re- 
ally efficient in removing sources with light curves that are suf- 
ficiently different from a supernova light curve. For example, an 
active galactic nuclei (AGN) with a light curve very similar to a 
supernova will of course not be removed by this cut. Using this 
rejection scheme allows us to perform an automatic cut that will 
only reject extreme outliers (0.1%) in the supernova population. 
Additional manual inspection is needed to safeguard against pos- 
sible spurious candidates that have light curves similar enough 
to SNe to end up with likelihoods higher than the 99.9% thresh- 
old. This process is described, along with our final supernova 
candidates, in Section|5] 

Some of the problems related to supernovae with non- 
standard light curves may be avoided by using the technique of 
"fuzzy" templates, pioneered by Ro dney & Tonrvl {2009). This 
method enables the templates to cover more of the parameter 
space by giving the templates themselves an uncertainty ("fuzzi- 
ness"). In this work we did not use this technique but it can be 
implemented into our code, and we will include it as an option 
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in future versions of it. Another way to improve typing of these 
SNe is to use a more extended set of templates and to update the 
existing ones. 

4.3. Systematic and statistical errors of typing 

We used the Monte-Carlo-simulated supernova light curves to 
find the expected misclassification ratios for the two main super- 
nova types at different redshifts. Typing of the simulated light 
curves was performed using the same code and assumptions as 
when typing the observed supernovae. This approach allows us 
to estimate the total errors on the supernova counts including 
both systematic errors, which arise from the assumed priors and 
possible systematic errors in the photometry; and statistical er- 
rors from the photometric statistical errors and the photometric 
redshift errors. 

We investigated the errors by studying the misclassification 
ratios, fioss and fgai„ for the simulated thermonuclear and core- 
collapse supernova light curves. The ratios are defined as 



fi„ss(TN) ^ 



flossiCC) 



Nt 



NappiTN) 
NappiCC) 



fgainiTN) ^ 



Nr 



NappiTN) 

Ntn^cc 

Nappiccy 



(17) 



where Ntn^cc and Ncc^tn is the number of SNe that were in- 
correctly classified (the true type on the left) and Ngpp is the ap- 
parent number of SNe of the given type. For example, fiossiTN) 
gives the fraction of thermonuclear SNe incorrectly classified as 
CC SNe. Note that using the ratios as measures of misclassifica- 
tion errors is valid if the underlying core-collapse to thermonu- 
clear SN ratio in the survey is the same as the one we assume in 
Sectionl4TT] 

The different subtypes have quite different misclassifica- 
tion ratios. Because we assumed the SN frequency for the 
different subtypes to be equal, we also computed weighted 
misclassification ratios, fhssATN), fgain.w(TN), fioss,w(CC), 
fgaiiuwiCC), where the number of misclassifications per subtype 
(e.g. Njip^tn) are weighted by the relative frequency of that sub- 
type. The subtype fractions are given in Table |2] The fioss.w and 
fgain.w ratios Were adopted as the final measure of misclassifica- 
tion errors in our survey. 

Figure [3] shows how fgain,w(TN/CC) and fioss,w{TN/CC) de- 
crease when a limit on the type probabilities, P/„„, for a source to 
be included in the analysis is introduced. The figure also shows 
that the use of a limiting probability will decrease the survey ef- 
ficiency. Owing to the low-number statistics nature of our survey 
we elected not to use a limiting probability as a rejection crite- 
ria, but it should be noted that the misclassification errors can be 
somewhat decreased by only including objects with high proba- 
bilities. The ratios and total number of simulated supernovae for 
the different redshift bins can also be found in Table |3] 

There is a caveat with using the misclassification results from 
the simulated sample. The errors are only really reUable when 
the observed supernova sample has light curves that are similar 
to the template light curves on average. If most of the observed 
supernovae have light curves different from the templates, the 
errors estimated by this method will not be representative. Figure 
ID shows the distribution of Bayesian evidence for most likely 
type, X, for the simulated SN sample together with the observed 
supernovae (see Section |5]l. Comparing the distribution in X of 
the simulated SNe to the observed ones indicate that the samples 
seems to be fairly different. 



A K-S test applied to the two samples is not conclusive, we 
cannot rule out that the two populations are the same with sig- 
nificance (the K-S p-value is 0.11). Interpreting this result is not 
trivial, there is a number of assumptions that go into simulating 
the SNe (the choice of specific templates and basically all of the 
priors) that can make the distribution different from what you 
would expect for a real sample. It is certainly possible that a dis- 
covered SNe may have a light curve that is simply different from 
any of the templates we use, which will cause a mismatch like 
the one we see in the X. distributions. This is also discussed in 
Section l4~2l where we provide some suggestions on how to solve 
the problem. 

The assumption that the the supernovae will be evenly dis- 
tributed into the subtypes (see Section 1431 1 will influence the X. 
distribution and cause it to be slightly different from the observed 
distribution. In this particular case we cannot use the a priori 
subtype fractions to weight the results (as used in Section 14.3b 
because the numbers of observed SNe in the different types are 
too low for subtype K-S calculations. 

There is consequently no reason to expect that the the sam- 
ples would be perfectly matched (giving p-values close to 1). 
Nevertheless, we cannot rule out that the resulting p-values indi- 
cate that the assumption of the simulated light curves are repre- 
sentative of the sample of real supernovae (at a significance level 
of 10%). Note that the efficient number of data points for the two 
types are quite low (~ 7 for the thermonuclear and ~ 9 for the 
core-collapse SNe, see Section |5]l. A better accuracy for the p- 
values requires more real supernovae to compare with, which 
would also enable us to look at the subtype statistics. 

4.4. Redshift determination with the typing code 

By fixing the type to the most likely template found in the full 
Bayesian fitting run, we tried to find the most likely redshift of 
the supernovae (see Section lJSj for a more detailed description). 
The prior on redshift is based on the probability distribution for 
photometric redshift of the host galaxies. By looking at the most 
likely redshifts obtained from the typing of the simulated su- 
pernovae, we investigated whether the prior information added 
through the supernova light curve changes the redshift accuracy. 

For the SVISS host galaxies the normalised photometric 
redshift scatter is cr, = rms((z,„(e - Zofo)/(l + Ztrue)) = 0.06, 
the simulated SN redshifts are scattered according to this (see 
Section 14.1b . When comparing the most likely redshifts {zfu) 
from the typing code with the true ones {ztme as defined in 
Section l4m i. we found a normalised redshift scatter of 0.067 for 
the thermonuclear supernovae and 0.072 for the core-collapse 
with a negligible offset for both types. Overall, a small number 
of objects (<0. 1 %) were assigned a catastrophic redshift, defined 
as objects with (|z/,r-Zm(cl/(l +Ztrue) > 0.2; these are all misclas- 
sified SNe. The resulting redshift accuracy is thus very similar 
to the input simulated error, the inclusion of prior information 
in the form of the supernova data did not improve the redshift 
determination. On the other hand, excluding the misclassified 
supernovae, no additional errors to the redshift estimates seem 
to have been introduced in the typing code. 

In Table [3] we also present the redshift scatter for the simu- 
lated supernovae in the redshift bins used to study misclassifica- 
tion. Note that the scatter values in the table are not normalised 
(i.e. dz - stdev(z/n(f - Zobs))), as opposed the the cr- discussed in 
the previous section. The overall redshift accuracy is consistent 
with the simulated input accuracy. The 0.25 < z < 0.5 redshift 
bin has a somewhat increased scatter, this is because of a higher 
contribution of catastrophic redshifts (or, equivalently, misclas- 
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Fig. 3. Effect on the misclassification ratios of applying a rejection criterion on the type probability. The left/right panel shows the 
weighted misclassification ratios {fioss,w and fgain.w) for the thermonuclear/core-collapse SNe along with the fraction of rejected 
candidates (Nrej/Napp) for a given limiting type probability {Pum)- This information can be used to select an optimal P/,„, to avoid 
misclassifications without rejecting too many candidates. For details see the text and Table|3](this figure is available in colour in the 
electronic version of the article). 

Table 3. Simulation results 



z 




fgain.w 


floss,w 


fga'm 


floss 


dz 


N 


Thermonuclear Supernovae 


z < 0.25 




0.019 


0.025 


0.039 


0.046 


0.091 


540 


0.25 < z 


<0.50 


0.098 


0.023 


0.096 


0.043 


0.12 


2894 


0.50 < z 


<0.75 


0.096 


0.032 


0.10 


0.052 


0.084 


3102 


0.75 < z 


< 1.00 


0.057 


0.017 


0.078 


0.022 


0.074 


1759 


1.00 < z 


< 1.25 


0.038 


0.0045 


0.061 


0.0056 


0.095 


358 


1.25 <z 


< 1.50 


0.26 


0.0 


0.50 


0.00 


0.082 


20 


z< 1.5 




0.082 


0.024 


0.042 


0.071 


0.098 


8673 


Core-collapse Supernovae 


z < 0.25 




0.0079 


0.0060 


0.015 


0.013 


0.078 


1676 


0.25 < z 


< 0.50 


0.020 


0.085 


0.038 


0.084 


0.12 


3319 


0.50 < z 


< 0.75 


0.041 


0.12 


0.067 


0.13 


0.083 


2428 


0.75 < z 


< 1.00 


0.026 


0.085 


0.032 


0.12 


0.094 


1172 


1.00 < z 


< 1.25 


0.0034 


0.029 


0.0042 


0.047 


0.12 


472 


1.25 < z 


< 1.50 


0.0 


0.014 


0.00 


0.027 


0.16 


367 


z< 1.5 




0.022 


0.075 


0.084 


0.037 


0.11 


9434 



Notes, fgain.w/fioss.w SlTS the misclassification ratios weighted by subtype frequency and fgami floss are the unweighted misclassification ratios, dz is 
the redshift variance (based on the fitted redshift from the typing code compared to the true redshift of the supernova) for the specific bin and N„,,p 
is the apparent number of simulated supernovae of the given type in the bin. 



sified SNe) in this bin. The actual percentage of outliers in this 
bin is still quite low, 3% for the thermonuclear and 4% for the 
core-collapse SNe. 

4.5. SDSSII supernovae 

We also tested our typing code on a small sample of spec- 
troscopically confirmed supernovae from the SDSS supernova 
survey (Frieman et al. 2008). The sample contains 55 la super- 
novae at z^O.001-0 2 and 32 IIP supernovae at z^O.OO 1-0.2 
dFrieman et al.l2008HKessler et al.l2009ll2010btlD' Andrea et alJ 



I20T0I) . The observations were obtained with a cadence (1-10 
days) and total survey length (~90 days). To be able to compare 
the typing results from this sample with the observations and 
simulations for the SVISS, we re-sampled the SDSS supernova 
light curves to a cadence of ~20 days (yielding 4-5 epochs). The 
resulting light curves (at z ~0.1) therefore sample approximately 
the same rest-frame epochs as the light curves from SVISS (with 
a mean redshift of 0.57). We used the SDSS g and r filters, which 
at z ~Q- 1 target similar rest-frame colours as the VIMOS R and 
/ filters at z ~ 0.5. Using only these epochs and filters will of 
course render the typing less optimal, but will render the re suits 
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Fig. 4. Main panel: fit quality as measured by the absolute Bayesian likelihood for the best-fitting template iX.max), versus P, the 
probability of a given source being of the simulated type, for the ~ 18000 Monte-Carlo-simulated supernovae. The upper panel 
shows the P distribution and the right panel the X-max distribution. For clarity, large (blue) dots and solid (blue) lines marks the 
correctly typed SNe (with P > 0.5); small (red) dots and dashed (red) lines the incorrectly typed SNe. The cross and star symbols in 
the main panel show X.max and P{TN)/P{CC), respectively, for the observed supernovae (see SectionlJJ. Three of the thermonuclear 
SN candidates have -CmoA < 1 x 10"'"^, they are marked by a caret lower limit symbol at X.max < 1 x 10"'"* (note that they have the 
same P(TN) values, thus the symbols are overlapping) (this figure is available in colour in the electronic version of the article). 



comparable to the SVISS typing because the same rest-frame 
light curve information is used. Note that a maximum of five 
epochs can be fitted using this sampling, compared to the maxi- 
mum of seven used in the SVISS observations and simulations. 

The results of testing performed on the sample of SDSS 
supernovae do not show any significant differences with the 
Monte-Carlo-simulated sample in terms of misclassification. 
The total misclassification percentage for a sample of 55 la su- 
pernovae (i.e., la supernovae typed as any of the core-collapse 
subtypes) is 7.3 %. The total misclassification percentage for a 



sample of 32 IIP supernovae is 9.3 %. These ratios are of the 
same order of magnitude as the results obtained for the simu- 
lated sample, although slightly higher The simulated la SNe are 
misclassified in ~3% of the cases in the relevant redshift bin, the 
corresponding percentage for IIP SNe is ~5%. The difference 
could be caused by the somewhat lower number of data-points 
available, on average, for the SDSSII SNe with resampled light 
curves. We conclude that the simulations provide valid estimates 
on the misclassification ratios for a real sample of supernova 
light curves. 
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Fig. 5. R-band image cutouts, approximately 20x20"large of the supernova candidates SVISS-SN43, SVISS-SN161 , SVISS-SNl 15 
and SVISS-SNl 16. The leftmost panels show the reference image, middle panels show the peak brightness (observed) epoch and 
the rightmost panels show the subtracted image at peak brightness. The (red) circle marks the location of the supernova as detected 
in the subtracted frame. The most likely redshifts (from the typing code) for these SNe are (starting at the top) 0.43, 0.50, 0.40 and 
0.55, respectively. 
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Fig. 6. R-band image cutouts for the supernova candidates SVISS-SN309, SVISS-SN402, SVISS-SN135 and SVISS-SN14. See 
Figure |5] for descriptions on the different panels and marks. The most likely redshifts (from the typing code) for these SNe are 
(starting at the top) 0.47, 0.22, 0.98 and 0.36, respectively. 
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Fig. 7. R-band image cutouts for the supernova candidates SVISS-SN51, SVISS-SN54, SVISS-SN261 and SVISS-SN55. See 
Figure |5] for descriptions on the different panels and marks. The most likely redshifts (from the typing code) for these SNe are 
(starting at the top) 0.51, 0.77, 0.37 and 0.83, respectively. 
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Fig. 8. R-band image cutouts for the supernova candidates SVISS-SN31, SVISS-SN56, SVISS-SN357 and SVISS-SN24. See 
Figure |5] for descriptions on the different panels and marks. The most likely redshifts (from the typing code) for these SNe are 
(starting at the top) 0.12, 0.57, 1.4 and 0.81, respectively. Note that the host galaxy for SN357 is blended with a foreground galaxy, 
the selected host is the lowermost of the two galaxies. 
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Fig. 9. Observed and fitted light cui-ves for the supernova candidates SVISS-SN43, SVISS-SN161, SVISS-SN115, and SVISS- 
SNl 16. The left-hand panels show the observed R (solid hexagons) and / (squares) light curves along with the best-fit light curve 
of the most likely supernova subtype, solid (black) for R and dashed (blue) for I. The error bars given for the observations are based 
on the photometric accuracy simulations described in Section 12.3.21 If the source is non-detected (i.e., has an estimated magnitude 
error of more than 1), a magnitude lower limit is given. The right-hand panels show the R - I colour evolution. For the R - I plot, 
the limit symbols indicate that the source was only detected in one of the bands and a lower or upper limit is given (this figure is 
available in colour in the electronic version of the article). 
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Fig, 10. Observed and fitted light cui-ves for the supernova candidates SVISS-SN309, SVISS-SN402, SVISS-SN135, and SVISS- 
SN14. See Figure|9]for descriptions on the different panels and marks. 
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Fig. 11. Observed and fitted light curves for the supernova candidates SVISS-SN51, SVISS-SN54, SVISS-SN261, and SVISS- 
SN55. See Figure|5]for descriptions on the different panels and marks. 
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Fig. 12. Observed and fitted light curves for the supernova candidates 
SN24. See Figure|5]for descriptions on the different panels and marks. 
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Table 4. Supernovae in the SVISS 



SVISS ID 


Type (Subtype) 


PAP.r 


log£ 


z 


Ms.peak 


te (days)" 


Ext./£(s - vy 


SVISS-SN43 


TN(Ia-normal) 


0.886(0.733) 


4.7 


0.43 


-19.88 


140 


-/o.o 


SVISS-SN161 


TN(Ia-faint) 


1.00(0.910) 


-19.1 


0.50 


-18.28 


-47 


-/o.o 


SVISS-SN115 


TN(Ia-91t) 


1.00(1.00) 


-17.8 


0.40 


-19.04 


1 


-/o.o 


SVISS-SN116 


TN(Ia-faint) 


1.00(1.00) 


1.1 


0.55 


-18.98 


-6 


-/o.o 


SVISS-SN309 


TN(Ia-faint) 


1.00(1.00) 


-14.8 


0.47 


-17.96 


-77 


-/o.o 


SVISS-SN402 


TN(Ia-faint) 


1.00(0.929) 


6.1 


0.22 


-18.66 


-127 


CZ/0.4 


SVISS-SN135 


TN(Ia-nomial) 


0.950(0.580) 


9.6 


0.98 


-18.9 


-25 


C3/0.2 


SVISS-SN14 


CC (Iln) 


0.623(0.551) 


11.2 


0.36" 


-18.64 


133 


-/O.O 


SVISS-SN51 


CC (IIP) 


1.00 (0.980) 


9.0 


0.51 


-17.56 


25 


C2/0.6 


SVISS-SN54 


CC (Iln) 


0.812 (0.733) 


9.8 


0.77 


-19.04 


72 


C2/0.2 


SVISS-SN261 


CC (Ibc-normal) 


0.734 (0.364) 


10.1 


0.37" 


-17.48 


-110 


C2/0.4 


SVISS-SN55 


CC (IIP) 


0.995 (0.854) 


3.8 


0.83 


-17.91 


-17 


-/O.O 


SVISS-SN31 


CC (IIL) 


0.999 (0.795) 


4.4 


0.12" 


-16.52 


64 


-/O.O 


SVISS-SN56 


CC (Ibc-bright) 


0.930 (0.912) 


0.037 


0.57 


-19.27 


-74 


C 1/0.6 


SVISS-SN357 


CC (Iln) 


1.00(1.00) 


6.0 


1.4 


-19.08 


-91 


-/O.O 


SVISS-SN24 


CC (Ibc-bright) 


0.643 (0.643) 


7.2 


0.81 


-19.56 


25 


C2/0.5 



Notes. P, refers to the co-added probability for the best-fitting main type (TN/CC), refers to the subtype probability from the typing code. 
**' This is the time since explosion in the observers frame. '''' Extinction models are as follows: C3, Cardelli with Rv=3.l; C2, Cardelli with Rv=2. 1 ; 
CZ, Calzetti. For the SNe with a most likely extinction of zero, no best-fit model is given. *"' "Hostless" SNe, for these a flat redshift prior was 
used, see text for details. 



5. Results 

After applying the photometric rejection criterion described in 
Section l2.3.2l we ended up with 115 supernova candidates that 
are typed using our code. The subtype probabilities for the two 
main types are co-added, yielding P{TN) and PiCC), for each 
candidate; the candidate was assigned the main type with the 
highest probability. £.„,ax for the candidate is the maximum evi- 
dence among the subtypes belonging to the chosen main type. 
The evidence was then compared to the threshold evidence, 
■£99.9,7- (see Section l4~2l i. and candidates with lower £.,„ax than 
the threshold were rejected. The remaining number of candi- 
dates at this point was 54, the rejected objects were also man- 
ually checked to ensure that the technique was working. 

We then investigated each of the 54 possible supernovae 
manually, making an overall assessment - using both the im- 
ages, light curve and classification - of whether the source is 
likely a supernova or something else. This investigation was per- 
formed by six of us (JM, TD, GO, LMT, JS, and SM), and we 
individually and separately rated each candidate. The manual re- 
jection step allows us to remove spurious transient sources from 
imperfect image subtraction. These residuals can be present in 
both bands and in multiple epochs, which enabled them to get 
through the earlier rejection steps. During the manual inspec- 
tion of both images and light curves, the residuals can be dis- 
covered. The choice of having several independent inspectors 
minimises the risk of errors being made. The rejected residuals 
have a number of properties in common: they all have a bright 
host galaxy, in general they tend to move (by 0.1-1 pixels) from 
epoch to epoch, most of them have low Xmox and have erratic 
light curves combined with large errors, most of them also show 
a negative residual close to the source in at least one epoch/filter 
Approximately 20 candidates are considered to be spurious sub- 
traction residuals. A small number of candidates (< 5) is found 
very close to the edge of the images and are related to the the 
higher noise level at the edges. A reanalysis of the photomet- 
ric errors in these regions showed that the sources did not fulfil 
our photometric criteria, we consequently decided to reject these 
sources. This decision also means that the efficient field of view 
will be somewhat smaller, a trade-off we are willing to accept 



to make the detection efficiency of the survey constant over the 
full field. After this rejection step we were left with 31 transient 
sources that we believe are real. The world coordinates, RI light 
curves and errors for the 31 transient sources are given in the 
appendix (available in the electronic version of the article). 

At this point we also rejected possible AGN that contaminate 
our sample. To do this we used subtracted frames with a two- 
year difference in time (the reference epoch from August 2003 
and the control epoch obtained in November 2005 to January 
2006). If a candidate showed variability over this time span, it 
is very unlikely to be a supernova, and we rejected it (with the 
exception of one candidate, discussed in Section 15.2b . Of the 
31 candidates 15 were rejected because of this, but note that 
some of these are likely not AGN but rather some other non- 
SN transient object (at least one of them has a light curve and 
colour consistent with a variable star). It should be noted that 
the AGN rejection scheme will not allow us to get rid of AGN 
that show no variation over the two-year baseline and have a 
SN-like light curve during the search period; however, the num- 
ber of AGN fulfilling this criterion is estimated to be small (also 
see Section l53] l. 

The final sample of 16 supernovae is presented in Table |4] 
cutout images for each SN are available in Figures |5]through|8] 
Figures l9]through[T2]show the observed and best-fit light curves 
for all of them. The best-fit light curves are here obtained as de- 
scribed in Section I3T2I it should be noted that the typing method 
is not based on finding the best-fitting (in the sense) param- 
eter set but rather on finding the type with the closest matching 
parameter space given the priors. Nevertheless, the fitted light 
curves provide a quick estimate on how good the evidence for 
the most probable type is, a fitted light curve far from the ob- 
served curve usually also means that the evidence is low. The 
probability (P) might be close to one even for low £.„,ax, this re- 
flects that the chosen type might be improbable in the absolute 
sense but that the type is the best among the available alterna- 
tives. 

The redshift distribution, based on the fitted redshifts from 
the typing, of the supernovae is shown in Figure [13] The joint 
distribution (the filled bars in the plot) has a mean of 0.58. The 
thermonuclear supernovae are on average found at lower red- 
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shifts compared to those of the core-collapse (with a mean of 
0.51 compared to 0.64). The requirement that the supernovae 
should be detected in both filters makes this expected because 
thermonuclear SNe become very faint at redshifts > 0.8 in the 
7? -band owing to k-corrections. 

5.1. Thermonuclear supernovae 

Seven of the 16 supernovae in the survey are found to have a 
thermonuclear origin. Based on the simulation results, approx- 
imately 0.7 (10% of the total number) could be misclassified 
core-collapse supernovae. Inspecting the light curves and look- 
ing at the fit quality {P(TN) and X.best) we can see that one of the 
seven supernovae (SVISS-SN43) has a somewhat insecure typ- 
ing. The probability, P{TN), is well below 1 .0, which means that 
it has a significant chance of being a core-collapse supernova. 

Three of the supernovae (SVISS-SN161, SVISS-SNl 15, and 
SVISS-SN309) have low evidences {X.best)- They all have prob- 
abilities close to 1.0, however, indicating that they, while not be- 
ing perfectly matched by the available templates and priors, are 
optimally fit by the thermonuclear templates. The addition of a 
more general stretch prior as described in Section |3]would most 
probably allow the templates to match the observed light curves 
better. 



5.2. Core-collapse supernovae 

Nine of the supernovae are found to be best matched by a core- 
collapse template light curve. The simulation results indicate 
that none of these are likely to be misclassified thermonuclear 
supernovae (2.4% of the total number). A look at the actual fit 
qualities for the supernovae show that three of them (SVISS- 
SN14, SVISS-SN261, and SVISS-SN24) have probability val- 
ues lower than 0.8 and accordingly have a definite probability of 
being misclassified thermonuclear SNe. It is worth noting that 
for two of them the most likely subtype is Ibc, the CC subtypes 
with light curves most similar to the thermonuclear subtypes. 



g3 
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Fig. 13. Redshift distribution of the supernovae detected in the 
SVISS ELAIS-Sl field using a z bin size of 0.25. The redshifts 
used are the resulting best-fit redshifts from the supernova typing 
code. The grey, filled bars show the total supernova counts, the 
blue (dotted) lines show the distribution of thermonuclear SNe 
and the red (dashed) lines the distribution of core-collapse SNe 
(this figure is available in colour in the electronic version of the 
article). 



They are also all quite faint and detected only in 2-3 epochs. 
However, we have no reason to mistrust the likelihoods com- 
puted for these SNe and thus conclude that they are likely of 
core-collapse origin. 

One of the CC supernovae, SVISS-SN56, has a low evi- 
dence (also visible in the light curve fit, see Figure [12), indi- 
cating that none of the supernova templates fitted the data well. 
Including more core-collapse supernova templates, improving 
the ones used or adding a colour prior could make it possible 
to get a better typing result for this supernova. 

As discussed in Section lZ4l for three CC SNe no photomet- 
ric redshift prior was used. In the case of SVISS-SN14 this is 
because the most likely redshift (0.36) for the found subtype, 
Iln, is significantly different from the photometric redshift of the 
host galaxy (0.76). This type and redshift is found to be the best 
fit even when the photometric redshift is used as a prior The 
reason for this mismatch could be either that the host redshift 
is incorrect (approximately 10% of the photometric redshifts are 
so called catastrophic redshifts, off by more than 0.3, again see 
Section 12.41 1. or that there is another very faint and undetected 
galaxy that is the real host of the SN. A third possibility is that 
the supernova is of a peculiar type that is not well-modelled by 
the template light curves. 

For SVISS-SN261, the closest galaxy is quite far away in 
angular distance and it has a redshift of 1.64. Using this as a 
redshift prior results in the SNe being placed at a redshift of 0.37 
as a Ibc-bright subtype. There is just no possible SN type that 
can give the observed light curve at such a high redshift. Also, 
using a flat prior on the redshift results in the same type and 
redshift. We believe that the true host of this SN is undetected, 
it could possibly be a low-surface brightness, but star-forming 
galaxy at intermediate redshift. 

The third of the three "hostless" SNe is a complicated case. 
SVISS-SN31 is located very close to a very faint galaxy with a 
photometric redshift of 2.97. At that redshift all of the supernova 
types under consideration would be undetectable (owing to dis- 
tance dimming and k-corrections). Using a flat prior on the red- 
shift results in a type IIL SN at z = 0.12. Other nearby galaxies 
(in angular distance) have been considered as hosts, but none of 
them has a redshift lower than ~0.6. The second-closest galaxy 
has a photometric redshift of 1.04 (highly uncertain, the 68% 
confidence interval for this measurement is 0.96-1.58), and if 
this is used as a prior for the supernova typing the resulting type 
is a Iln SN at z = 1.92. However, the evidence for this type is 
lower than the IIL solution found with a flat prior Redshift priors 
from the other nearby galaxies result in even worse evidence for 
the most likely types. It should also be noted that the main type 
(CC) of the supernova is the same for all the tested redshift pri- 
ors. Based on the evidence comparisons, host galaxy properties, 
and angular distance we conclude that the typing without a red- 
shift prior provides the most believable solution. The mismatch 
with the host redshift could be explained by, as for SVISS-SN14 
above, the host galaxy redshift being catastrophically wrong, a 
faint undetected host galaxy, or a peculiar SN type. 

5.3. A peculiar SN candidate 

In Section|5]we described how the AGN that likely contaminate 
the sample are dealt with. One of the variable sources showing 
variability over a two year period, SVISS-T167 has an intrigu- 
ing light curve that could possibly come from a supernova with 
a very extended light curve. However, the supernova would in 
that case have to be of a peculiar type; our typing code was un- 
able to find the most likely type (i.e. Jlhest < -€99.9, as defined 
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Fig. 15. /-band image cutouts, approximately 20 x 20"large of the transient of unknown type, SVISS-T167. The leftmost panel 
shows the reference image, middle panel shows the peak brightness (observed) epoch and the rightmost panel shows the subtracted 
image at peak brightness. The (red) circle marks the location of the transient as detected in the subtracted frame. 



PISNo@z=l.O 



Epoch 

Fig. 14. /-band light curve of SVISS-T167. Also shown is a B- 
band light curve for a pair instability SNe shifted to z = 1 (line). 
Details of how the comparison light curve is constructed can be 
found in the text (this figure is available in colour in the elec- 
tronic version of the article). 



in Section I4.2l i. The light curve of this candidate is shown in 
Figure[T4]and an image cutout is shown in Figure[T5] For it to be 
even remotely similar to a normal supernova light curve, it needs 
to be situated at high redshift where cosmological time dilation 
will stretch the light curve by a factor of 2 (or more). 

Comparing this to the light curves of bright peculiar super- 
novae we find that the light curve is marginally consistent with 
that of a pair instability (PI) supernovae at high redshift, some- 
thing that has also been suggested by Gal- Yam et al. (2009) for 
SN 2007bi and by Woosley et al.1 (12007) for SN 2006gy. PI SNe 
have slowly evolving light curves, sometimes with a secondary 
peak; the slow light curve decline of our candidate indicates that 
it could be similar to such SNe. The PI light curve in Figure fl4l 
was derived by taking a model rest frame B band light curves for 
a 110 Mq progenitor from the supplemental information avail- 
able for IWoosIev et al.l (l2007h and shifting it to z = 1 .0, ap- 
plying time dilation and distance modulus. We did not attempt 
any detailed k-correction, but it should be noted that the correc- 
tion term should be fairly small because the rest frame B central 
wavelength at z - 1 .0 is close to the central wavelength of the 
/ filter We also compared the light curve to that of SN 2010gx 



(iPastorello et al.ll2010h . an extremely luminous type Ic-like su- 
pernova. We found that the very extended light curve and colour 
evolution of our transient is impossible to match with the 2010gx 
light curve, even when shifting it to higher redshift. 

There is no strong evidence for our candidate being a super- 
nova of this type, it is still a distinct possibility that the source 
is actually an AGN because the position of the source is right in 
the centre of its host galaxy. The photometric redshift for the host 
galaxy is not constrained, which might point towards it being an 
AGN (our photometric redshift code does not include AGN SED 
templates), but could also mean that it is a high-redshift galaxy 
(z > I) where the lack of near infrared data renders the pho- 
tometric redshift technique unreliable. We excluded the source 
from the SN sample presented in Table |4] To discover the true 
identity of the object, follow-up observations of its host galaxy 
are needed. 

6. Summary and discussion 

We have presented light curves and typing of the 16 supernovae 
discovered by SVISS. Using photometric redshifts estimated 
from host galaxy photometry, we typed the SNe into thermonu- 
clear or core-collapse events. The typing uses spectral and light 
curve templates for nine different subtypes and is based on a 
Bayesian fitting method. From studying Monte-Carlo-simulated 
SN light curves and a small sample of SDSS SN light curves, we 
conclude that the typing code yields reliable results considering 
that we only have photometric redshifts with accuracies in the 
order of 0.1, only one colour, and a limited number of epochs. 
We find that probably not more than 5-10% of the supernovae 
are misclassified. The misclassification ratios are not symmetric 
with respect to the two main types, a CC SN is more likely to be 
misclassified into a TN SN than the other way around. 

Comparing our typing quality with that of other typing 
codes is not straight-forward because the input data and tests of 
the codes vary considerably. iKessler et al.l (1201 Oah try to rem- 
edy this by offering a test sample of SN light curves and a 
framework in which different typing codes can be compared. 
Currently, our code is written specifically to deal with the data 
from SVISS, trying the code out on the test sample is thus 
outside the scope of this paper Nevertheless, we can compare 
the typing quality estimated from our Monte-Carlo-simulated 
light cu rves with the typing qualities reported in Kessler et al] 
(1201 Oal) . The la efficiency, which in our paper is equal to 
{NappiTN) - Ntn^cc)/N,o,{TN), and lapureness, {NappiJN) - 
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NTN^cc)/(Napp(TN) - Ntn^cc + Ncc^tn) in ou r notation, are 
two of the concepts used to describe the quality in'Kessl er et alJ 
(1201 Oah . Our overall efficiency and pureness values for the ther- 
monuclear SNe, 0.97 and 0.75, respectively, are comparable to 
the results obtained by the other typing codes, but with the caveat 
that we did not use the same test sample of supernovae as these 
authors. 

Out of the 16 SNe, 9^q2 are core-collapse supernovae, and 
are thermonuclear; only the systematic typing errors de- 
rived from the misclassification ratios are given here. The mean 
redshift for the entire sample is 0.58 with the core-collapse SNe 
being at slightly higher redshifts than the thermonuclear ones. 
In a supernova survey of this type a large number of spurious 
sources will be detected and methods must be devised to safely 
reject interlopers without compromising the SN detection effi- 
ciency. We used a combination of automatic (SExtractor, signal- 
to-noise rejection, SN typing rejection) and manual (by-eye in- 
spection of image frames by multiple independent observers, 
by-eye inspection of light curves by multiple independent ob- 
servers) rejection methods to ensure that the final SN sample is 
as pure as possible. 

In the specific case of AGN interlopers we find that 
the use of an additional observing epoch displaced from the 
search/reference epochs by at least one year makes it possible to 
limit the amount of contamination further We also find that most 
of the AGN will also be rejected in the supernova typing step, be- 
cause the fitting code is simply not able to fit the SN templates to 
their colour and evolution. At least one of the rejected candidates 
could possibly be a supernova of an extremely rare type. This 
particular source was still detectable, albeit barely, one year after 
our search period ended. The extremely e xtended light curve is 
comparable to that of pair instabiUty SNe (IWooslev et al.ll2007h 
at high redshifts. 

The SVISS results show that SN imaging searches combined 
with photometric typing can be used to find and characterise SNe 
at high redshifts. While the basic strategy used in this project 
does work, there are a few improvements that could be con- 
sidered. One of the main uncertainties in the SN typing is the 
use of photometric redshifts. Securing more precise redshifts for 
the host galaxies would make the typing considerably more ac- 
curate. These redshifts could be in the form of well-calibrated 
photometric or spectroscopic redshifts. Future pencil beam sur- 
veys targeting high-redshift SNe should thus be aimed at regions 
with good spectral and multiwavelength coverage. We found the 
cadence of the search to be more or less optimal, given the 
constraint on total observational time, but with more frequent 
observations it should be possible to det ect the UV breakou t 
phase of core collapse SNe as well (e.g., Tominaga et al.ll20TlT) . 
Observations in more than one filter is mandatory to be able to 
do photometric typing, but we also found it to be quite important 
for rejecting spurious objects in the detection step. Adding more 
filters will break some of the degeneracies in the typing, but the 
typing accuracy is also limited by the quality of the templates. 

A full analysis of supernova and star-formation rates based 
on our observations will be presented in a forthcoming paper 
(Melinder et al., in prep). The work presented in this paper will 
help us estimate the uncertainties of the found rates. To calcu- 
late proper rates, information about the detection efficiency and 
the enforced photometric limits has to be taken into account. 
Nevertheless, a quick look at the raw numbers of discovered 
SNe of the two types compared to Monte-Carlo-simulations of 
detectable SNe using supernova rates based on local o bserva- 
tions and star-formation rates (Dahle n & Franssonll 1999b shows 
that our numbers are consistent with the expected rates. We can 



also compare the range and distribution in redshift of our su- 
pernovae to that of other surveys. Our sample of thermonuclear 
SNe a re not really uni que in this respect (see, e.g., Neill et al] 
120061: iPoznanski et alJ l2007b: D ahlen et al. 2008), but it should 
be noted that the survey methodology we used is different from 
most other surveys where the requirement of spectra means 
that some of the fainter events may be missed. For the core- 
collapse SNe our samp le is similar to the observations presented 
in Dahl en et all ( l2004l) with 17 CC SNe at z < 0.7. We detect 9 
CC SNe out to z = 1.4, with most of them at z ~0.8, and this is 
therefore one of the deepest survey to date for CC SNe. 
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Appendix A: Observational data for the 31 transient 
objects 

This appendix contains a table with extended information for 
the 31 transient objects detected in the SVISS ELAIS-Sl field. 
The astrometry is accurate to within 0.4"(for details on the as- 
trometry see Mencia-Trinchant et al., in prep.). The photomet- 
rical errors are obtained as described in Section I2.3.2J For non- 
detections the magnitude is given as a lower limit and the limit 
is also given as the error The limiting magnitudes given in the 
table are 1 cr limiting magnitudes derived from the photometric 
accuracy investigation. The host galaxy redshifts given in the ta- 
ble are photometric redshifts derived from the UBVRI galaxy 
photometry. Objects that were assumed to be hostless are given 
a redshift of N/A in the table. 
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Table A.l. Extended infomiation on the transient objects 



Id 



RA(deg.) Dec (dcg.) Host;„ 



)och 




Sm,i 


mi 


Siiii 


1 


24.881 


0.0811 


lA.lAl 


0.1405 


2 


25.416 


0.1149 


24.797 


0.1429 


3 


25.615 


0.1458 


25.47 


0.3558 


4 


26.427 


0.222 


25.518 


0.3195 


5 


26.649 


0.5808 


26.927 


0.5354 


6 


27.29 


0.421 


>26.894 


26.8935 


7 


27.36 


0.4706 


>26.719 


26.7194 


1 


>28.134 


28.1342 


>28.024 


28.0236 


2 


>27.94 


27.9396 


>27.746 


27.7464 


3 


25.115 


0.0595 


25.503 


0.1809 


4 


24.207 


0.02 


23.878 


0.0362 


5 


25.925 


0.1307 


25.383 


0.1815 


6 


27.213 


0.3706 


26.709 


0.4135 


7 


>28.084 


28.0837 


>27.722 


27.7219 


1 


>27.3 


27.3 


>27.211 


27.2112 


2 


22.07 


0.02 


22.235 


0.02 


3 


22.884 


0.02 


22.262 


0.02 


4 


24.644 


0.0676 


23.416 


0.02 


5 


24.965 


0.0839 


23.982 


0.0716 


6 


25.221 


0.1118 


24.703 


0.1808 


7 


25.36 


0.1697 


25.138 


0.3596 


1 


>28.657 


28.6567 


>27.598 


27.5984 


2 


22.895 


0.02 


22.941 


0.02 


3 


24.622 


0.1311 


23.45 


0.02 


4 


26.516 


0.2035 


25.441 


0.1349 


5 


26.335 


0.3015 


24.983 


0.1291 


6 


27.643 


0.5584 


25.911 


0.3007 


7 


26.733 


0.3231 


26.434 


0.3512 


1 


>28.657 


28.6567 


>27.598 


27.5984 


2 


>28.508 


28.5081 


>27.635 


27.635 


3 


>28.022 


28.0224 


>27.577 


27.5774 


4 


24.236 


0.02 


23.902 


0.0532 


5 


23.599 


0.02 


23.751 


0.0761 


6 


26.126 


0.1866 


25.278 


0.1809 


7 


27.559 


0.5768 


25.75 


0.2316 


1 


>27.3 


27.3 


>27.211 


27.2112 


2 


>26.776 


26.7764 


>27.543 


27.5435 


3 


>28.096 


28.0961 


>26.972 


26.9718 


4 


>27.703 


27.7034 


>27.298 


27.2979 


5 


>27.221 


27.2206 


>25.897 


25.8973 


6 


24.173 


0.02 


23.362 


0.02 


7 


23.747 


0.02 


23.063 


0.02 


1 


>27.253 


27.2533 


>27.877 


27.8767 


2 


>28.186 


28.1863 


>27.596 


27.596 


3 


25.344 


0.0522 


24.569 


0.0683 


4 


27.161 


0.3049 


25.469 


0.1887 


5 


27.815 


0.612 


26.092 


0.3977 


6 


>27.173 


27.1731 


26.736 


0.7313 


7 


28.259 


0.5435 


27.158 


0.5231 


1 


26.187 


0.1873 


25.843 


0.1986 


2 


26.291 


0.2177 


26.025 


0.2372 


3 


27.031 


0.3204 


26.979 


0.4307 


4 


27.948 


0.6347 


27.462 


0.5716 


5 


27.766 


0.5392 


27.096 


0.6211 


6 


28.03 


0.5945 


27.42 


0.5754 


7 


>28.084 


28.0837 


26.636 


0.3541 


1 


26.042 


0.2719 


24.929 


0.1633 


2 


26.663 


0.3074 


25.147 


0.1941 


3 


27.145 


0.4142 


25.47 


0.3558 


4 


26.9 


0.288 


25.967 


0.4471 



SVISS-SN43 7.9102451 -44.499318 0.46 



SVISS-SN161 7.9549322 -44.484425 0.54 



SVISS-SN115 8.1434494 -44.39665 



0.40 



SVISS-SN116 8.088403 



-44.47483 



0.57 



SVISS-SN309 8.0001281 -44.439905 0.67 



SVISS-SN402 8.0848318 -44.405671 0.34 



SVISS-SN135 7.8847856 -44.377629 0.90 



SVISS-SN14 8.0329574 -44.514378 N/A 



SVISS-SN51 7.9625753 -44.508399 0.53 
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Table A.l. Continued. 



Id 


RA (deg.) 


Dec (deg.) 


Host Zp 


Epoch 


niR 


dniR 


nil 


6mi 










5 


27.25 


OAllS 


25.418 


0.14 










6 


>28.19 


28.19 


25.954 


0.3464 










7 


>28.022 


28.0217 


26.682 


0.5789 


SVISS-SN54 


7.9628816 


-44.443032 


0.77 


1 


25.467 


0.16 


24.977 


0.1696 










2 


26.722 


0.3195 


25.563 


0.2678 










3 


>27.967 


27.967 


26.066 


0.6209 










4 


>28.386 


28.3863 


26.671 


0.6957 










5 


>28.275 


28.2753 


>28.148 


28.1477 










6 


>27.008 


27.008 


>26.489 


26.489 










7 


>28.022 


28.0217 


>26.719 


26.7194 


SVISS-SN261 


8.0816154 


-44.476649 


N/A 


1 


>28.134 


28.1342 


>28.024 


28.0236 










2 


>27.94 


27.9396 


>27.746 


27.7464 










3 


>28.285 


28.2846 


>27.675 


27.675 










4 


>28.171 


28.1708 


>27.643 


27.6432 










5 


>28.204 


28.2036 


>27.305 


27.3051 










6 


24.892 


0.0416 


23.8 


0.0426 










7 


26.361 


0.2134 


24.922 


0.1047 


SVISS-SN55 


7.8047583 


-44.499259 


0.59 


1 


>27.918 


27.9184 


>27.667 


27.6672 










2 


25.776 


0.1298 


25.566 


0.2291 










3 


25.878 


0.1712 


24.981 


0.1349 










4 


>27.418 


27.4183 


25.545 


0.3008 










5 


>28.065 


28.0654 


25.115 


0.2053 










6 


>28.133 


28.1333 


26.083 


0.3359 










7 


>27.753 


27.7532 


26.291 


0.4234 


SVISS-SN31 


8.0856037 


-44.462252 


N/A 


1 


24.303 


0.0459 


23.958 


0.0535 










2 


25.016 


0.0526 


24.313 


0.0634 










3 


25.671 


0.1325 


24.756 


0.1057 










4 


26.771 


0.248 


25.689 


0.151 










5 


27.031 


0.5771 


25.592 


0.2423 










6 


27.535 


0.5243 


26.24 


0.3787 










7 


>27.966 


27.9659 


26.308 


0.327 


SVISS-SN56 


8.0617464 


-44.438378 


0.51 


1 


>26.694 


26.694 


>26.359 


26.359 










2 


>26.776 


26.7764 


>26.333 


26.333 










3 


>28.096 


28.0961 


>25.939 


25.939 










4 


25.372 


0.1405 


24.312 


0.1262 










5 


25.343 


0.1372 


23.897 


0.0672 










6 


27.093 


0.5583 


24.884 


0.2127 










7 


>27.322 


27.3217 


>26.176 


26.1764 


SVISS-SN357 


7.9953784 


-44.396938 


1.36 


1 


>27.3 


27.3 


>27.211 


27.2112 










2 


>26.776 


26.7764 


>27.543 


27.5435 










3 


>28.096 


28.0961 


>26.972 


26.9718 










4 


>27.703 


27.7034 


>27.298 


27.2979 










5 


24.743 


0.0709 


25.001 


0.3677 










6 


24.922 


0.0813 


24.624 


0.1682 










7 


25.344 


0.1673 


24.729 


0.2138 


SVISS-SN24 


7.9023817 


-44.374341 


0.77 


1 


25.184 


0.138 


23.816 


0.0755 










2 


26.489 


0.2786 


24.369 


0.1448 










3 


>27.226 


27.2258 


26.055 


0.3293 










4 


>27.836 


27.8362 


26.635 


0.5763 










5 


>27.153 


27.1535 


>26.362 


26.3618 










6 


>27.451 


27.451 


>26.641 


26.6409 










7 


>27.503 


27.5032 


>26.398 


26.3982 


SVISS-T99 


8.0147036 


-44.388741 


1.50 


1 


26.438 


0.4123 


24.549 


0.1501 










2 


25.545 


0.3216 


24.117 


0.1073 










3 




25.096 i 


O /I TIC 


A 1 OAA 










4 


26.49 


0.329 


24.091 


0.1019 










5 


25.674 


0.216 


24.046 


0.077 










6 


25.861 


0.2149 


24.139 


0.1084 










7 


24.644 


0.1045 


23.611 


0.1222 


SVISS-T104 


7.9834014 


-44.414223 


1.52 


1 


24.255 


0.053 


23.965 


0.0922 










2 


24.567 


0.0901 


24.03 


0.0987 
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Id 


RA (deg.) 


Dec (deg.) 


Host Zp 


Epoch 


niR 




nil 


6mi 










3 


23.751 


0.02 


24.033 


0.0971 










4 


25.153 


0.1144 


24.243 


0.1186 










5 


26.532 


0.5599 


25.183 


0.468 










6 


27.25 


0.6157 


>26.62 


26.6198 










7 


>27.322 


27.3217 


>26.176 


26.1764 


SVISS-T167 


8.0149816 


-44.361213 


1.95 


1 


>27.3 


27.3 


>27.211 


27.2112 










2 


>16.116 


26.7764 


25.764 


0.2732 










3 


26.07 


0.2859 


25.337 


0.2675 










4 


25.587 


0.1695 


24.991 


0.2025 










5 


25.988 


0.3184 


25.02 


0.3775 










6 


25.193 


0.1084 


24.783 


0.1944 










7 


25.81 


0.2537 


25.215 


0.3938 


SVISS-T182 


8.1151531 


-44.431759 


0.52 


1 


>27.3 


27.3 


>27.211 


27.2112 










2 


>26.776 


26.7764 


>27.543 


27.5435 










3 


24.4 


0.2141 


24.431 


0.1453 










4 


23.198 


0.02 


23.078 


0.02 










5 


23.451 


0.02 


22.84 


0.02 










6 


23.332 


0.02 


23.095 


0.02 










7 


23.92 


0.02 


24.134 


0.1091 


SVISS-T23 


7.8220527 


-44.379293 


0.65 


1 


25.619 


0.254 


25.926 


0.5165 










2 


25.317 


0.1843 


25.601 


0.2831 










3 


25.48 


0.2386 


25.233 


0.2195 










4 


25.923 


0.2426 


26.226 


0.4586 










5 


25.259 


0.2237 


>26.362 


26.3618 










6 


26.155 


0.3105 


26.391 


0.446 










7 


26.25 


0.3229 


>26.398 


26.3982 


SVISS-T26 


7.9150698 


-44.334945 


1.23 


1 


>26.887 


26.887 


25.088 


0.234 










2 


24.828 


0.1365 


24.297 


0.1378 










3 


24.365 


0.1184 


23.805 


0.0795 










4 


24.65 


0.1249 


24.02 


0.0846 










5 


24.114 


0.02 


23.852 


0.0797 










6 


24.643 


0.111 


23.905 


0.3598 










7 


24.967 


0.1364 


24.367 


0.1622 


SVISS-T68 


7.8542558 


-44.400151 


0.34 


1 


>27.253 


27.2533 


>27.877 


27.8767 










2 


26.33 


0.193 


25.546 


0.1742 










3 


26.515 


0.3176 


25.981 


0.2312 










4 


26.989 


0.2703 


25.697 


0.2381 










5 


26.196 


0.1796 


25.744 


0.29 










6 


26.926 


0.8029 


>26.911 


26.9113 










7 


>28.472 


28.4721 


25.381 


0.159 


SVISS-T306 


7.9255526 


-44.517406 


0.4 


1 


23.315 


0.02 


23.44 


0.02 










2 


23.438 


0.02 


23.718 


0.0461 










3 


23.054 


0.02 


23.312 


0.02 










4 


23.795 


0.02 


24.143 


0.0783 










5 


23.294 


0.02 


23.281 


0.02 










6 


24.042 


0.02 


24.809 


0.153 










7 


25.272 


0.0876 


>26.719 


26.7194 


SVISS-T39 


7.9069841 


-44.529026 


1.8 


1 


22.849 


0.02 


22.691 


0.02 










2 


23.087 


0.02 


22.806 


0.02 










3 


23.366 


0.02 


23.124 


0.02 










4 


23.645 


0.02 


23.418 


0.02 










5 


24.422 


0.082 


23.565 


0.0596 










6 


24.115 


0.02 


23.516 


0.1144 










7 


22.971 


0.02 


22.845 


0.02 


C\7TCC T^CA 

aVlbJ>-l jU 


/.960/jzy 


-44.jl 161o 


yj.i 


1 


Zd./oZ 


A O 1 


ZJ.Z/Z 


A o 1 o/; 
U.Ziz6 










2 


26.444 


0.2649 


25.635 


0.2819 










3 


27.018 


0.3854 


>26.602 


26.6022 










4 


27.628 


0.4064 


>26.764 


26.764 










5 


>28.275 


28.2753 


>28.148 


28.1477 










6 


>28.19 


28.19 


>26.894 


26.8935 










7 


27.538 


0.5227 


26.654 


0.5722 
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Id 


RA (deg.) 


Dec (deg.) 


Host Zp 


Epoch 


niR 


dniR 


nil 


6mi 


SVISS-T6 


7.870627 


-44.513955 


1.64 


1 


24.408 


0.1208 


25.075 


0.2029 










2 


23.841 


0.02 


24.69 


0.1199 










3 


25.261 


0.1689 


25.071 


0.2404 










4 


23.709 


0.02 


24.144 


0.0849 










5 


24.466 


0.0946 


>26.476 


26.4764 










6 


>27.004 


27.0039 


24.518 


0.1963 










7 


24.187 


0.02 


24.602 


0.2067 


SVISS-T36 


7.8742713 


-44.409262 


N/A 


1 


23.593 


0.02 


22.525 


0.02 










2 


22.87 


0.02 


22.35 


0.02 










3 


23.355 


0.02 


22.776 


0.02 










4 


21.948 


0.02 


21.894 


0.02 










5 


23.726 


0.02 


21.124 


0.02 










6 


22.609 


0.02 


22.253 


0.02 










7 


23.133 


0.02 


22.636 


0.02 


SVISS-T59 


7.7726376 


-44.492906 


0.64 


1 


>27.364 


27.3641 


>26.725 


26.7249 










2 


27.552 


0.5635 


24.851 


0.1431 










3 


25.557 


0.2283 


24.719 


0.1742 










4 


>27.3 


27.2997 


24.528 


0.1315 










5 


25.371 


0.193 


24.148 


0.1936 










6 


24.868 


0.1861 


24.741 


0.2644 










7 


26.298 


0.3412 


24.968 


0.277 


SVISS-T106 


7.9025451 


-44.494797 


0.09 


1 


>27.364 


27.3641 


>26.725 


26.7249 










2 


>27.587 


27.5866 


>26.912 


26.912 
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